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Abstract. We investigate the behavior of various spectral invariants, particularly the de- 
terminant of the Laplacian, on a family of smooth Riemannian manifolds fl e which undergo 
conic degeneration; that is, which converge in a particular way to a manifold with a conical 
singularity. Our main result is an asymptotic formula for the determinant up to terms which 
kjj^ vanish as e goes to zero. The proof proceeds in two parts; we study the fine structure of the 

heat trace on the degenerating manifolds via a parametrix construction, and then use that 
fine structure to analyze the zeta function and determinant of the Laplacian. 



1. Introduction 

Let N be a. closed Riemannian manifold of dimension n — 1 with a fixed metric dy 2 . The 
infinite cone over N, which we call Cn, is the Riemannian manifold (0, oo) r x N y with the 
conic metric 

dr 2 + r 2 dy 2 . (1) 

Let Qq be a compact manifold of dimension n > 2 which has a single isolated conical 
singularity of cross-section A^ at a point P G fi . This means precisely that there is some 
5 > and some neighborhood of P which is isometric to (0, 5} r x N y , with the conic metric 
(|T]), where P corresponds to r = 0. Throughout, we assume without loss of generality that 
Q is exactly conic for r < 2. Further, let Z be a complete manifold, without boundary, 
which is exactly conic in a neighborhood of oo with cross-section N; that is, Z has the metric 
(00), but in a neighborhood of r = oo rather than r = 0. Again, without loss of generality, 
we assume that Z is exactly conic for r > ^. 

Using Z, we construct a degenerating family of smooth manifolds f2 £ that converge to Qq. 
For each e > 0, we cut Z at r = 1/e and then scale the metric by e 2 , obtaining the manifold 
e{Zfl{r < 1/e}}. We then replace the portion of f2o where r < 1 with e{Zn{r < 1/e}}. The 
metrics agree in a neighborhood of the boundary r = 1, so they may be glued together to 
obtain the smooth manifold Q e . This process is illustrated in Figured! Note that inside r = 1, 
Q t is isometric to eZ, while outside r = e, Q t is isometric to Qq. As e — > 0, the manifolds 
Q e converge in a natural way to Qq. We will investigate the behavior of the spectrum of 
the Laplacian, the trace of the heat kernel, and the determinant of the Laplacian on fl e as e 
approaches zero. 

1.1. Motivation. Families of smooth Riemannian manifolds which degenerate in some way 
to a manifold with conical singularities arise naturally in various settings in differential ge- 
ometry and PDE. One such setting is a program, initiated by Melrose and collaborators, 
to study elliptic differential equations on singular spaces through the methods of geometric 
microlocal analysis. The goal is to study the transition between smooth and singular geom- 
etry under the degenerating process, and in particular to examine the limiting behavior of 
various geometric and index-theoretic invariants. In this setting, conic degeneration was first 
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FIGURE 1. Conic degeneration: construction of 

introduced in the Ph.D. thesis of McDonald [McD] . who studied the behavior of the Schwarz 
kernel of the resolvent under a type of conic degeneration similar to the one we consider. 
The particular technique he used is known as analytic surgery. 

Along these lines, one of our ultimate objectives is a generalization of the Cheeger-Muller 
theorem to manifolds with conical singularities. The Cheeger-Muller theorem in the smooth 
setting concerns the analytic torsion, which is an invariant built out of an alternating sum of 
determinants of Laplacians, first defined in [RS] . The analytic torsion appears to be strongly 
dependent on the metric; however, the Cheeger-Muller theorem shows that it is in fact equal 
to a combinatorial invariant called the R-torsion. The theorem was first conjectured in [RSJ 
and then proven independently in [Chi J and [MiiJ. Analyzing the determinant of the scalar 
Laplacian under conic degeneration is a first step towards analyzing the limiting behavior 
of the analytic torsion. We hope to combine this analysis with an analysis of the R-torsion 
under conic degeneration to obtain a Cheeger-Muller theorem for manifolds with conical 
singularities. 

Analyzing the behavior of the determinant under conic degeneration may also have some 
bearing on the isospectral problem. In a well-known series of papers [OPS1, OPS2J, Osgood, 
Phillips, and Sarnak use the determinant of the Laplacian to prove that any set of isospectral 
metrics on a closed surface is compact. A key step in their approach is to show that the 
log of the determinant is a proper map on the moduli space of constant curvature metrics. 
Using similar techniques, they also prove isospectral compactness for planar domains in R 2 
with Dirichlet boundary conditions [ OP S3] . 

A natural question, asked by Khuri in |Khj . is whether these results generalize to sets 
of isospectral metrics on topological surfaces of genus p with n disks removed, again with 
Dirichlet boundary conditions, for np > 1. Khuri showed that the approach of |OPSlt lO~PS2, 
IOPS3] cannot work in this setting, because the log of the determinant is no longer a proper 
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map on the appropriate moduli space. In particular, Khuri constructed explicit families of 
constant curvature metrics which approach the boundary of moduli space but which have 
log determinants staying bounded from below. These families are constructed by taking 
a fixed flat metric on a surface of genus p with n conical singularities and then excising 
smaller and smaller disks around the conic points; see [KhJ for the details. This process 
is similar to our conic degeneration. We hope that understanding the precise behavior of 
the determinant under conic degeneration will help us gain a better understanding of the 
determinant on moduli space and provide new approaches to the isospectral compactness 
problem. We should also mention that Khuri's problem was recently partially solved by 
Kim, who proved isospectral compactness in the np > 1 case for flat metrics via a different 
approach [Kim] . 

1.2. Prior work. The particular form of conic degeneration we consider was first introduced 
in 2004 by Degeratu and Mazzeo |Ma2j . in slightly greater generality. Their goal was to 
analyze elliptic operators not just on manifolds with conical singularities but also on iterated 
cone-edge spaces. They studied the behavior of the eigenvalues of the Laplacian in this 
general setting, proving convergence of the spectrum under certain general assumptions 
[Ma2j . Our conic degeneration corresponds to the depth-1 case of their work, for which 
Rowlett proved the following spectral convergence [Ma2l IRowj : 

Theorem 1. [Degeratu- Mazzeo, Rowlett] Assume n > 5. Let X e ^ be the ith eigenvalue of 
the Laplacian on Q e , and let Ao,i be the ith eigenvalue of the Laplacian, with the Friedrichs 
extension, on Q . Then for each i, as e — > 0, 

A e ,j — » Ao,i- 

Rowlett also proves this result for n > 3 under some additional hypotheses |Row2j . The 
study of spectral convergence has been refined further by Anne and Takahashi in [ATJ. 

It turns out to be significantly harder to analyze the behavior of more complicated spectral 
invariants, such as the heat trace and determinant, under conic degeneration. For example, 
based on Theorem [TJ, one expects that the heat trace H n ^ (t) should converge as e goes to 
zero for any fixed positive t. However, this convergence is not easy to prove, because the 
convergence in Theorem [T] is not uniform in i. In general, heat trace convergence has only 
been studied under restrictive curvature assumptions; for example, if the sectional curvatures 
of Q e were uniformly bounded below (which is equivalent to the assumption that Z has non- 



negative sectional curvature), it would follow from |Ding . The assumption of non- negative 



sectional curvature is, however, quite strong, and we wish to avoid it. 

Mazzeo and Rowlett, in [MRj . study a closely related problem: that of degeneration of 
smoothly bounded planar domains in M 2 to domains with corners. This is illustrated in 
Figure [TT21 In this setting, with Dirichlet boundary conditions, spectral convergence holds, 
and as shown in |MRj . it is also true that for each positive t, 

TrH n '(t) ->TrH n °(t). (2) 

However, the fine structure of TrH Q,! (t) near t = e = is more complicated; §Zfj is not 
uniform in t. For each e > 0, the usual short-time heat asymptotics imply that there are a& j£ 
such that 

TrH^(t) = ao.et" 1 + a^r 1 ' 2 + a 2 , e + 0{t). 
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Figure 2. Polygonal degeneration of Mazzeo-Rowlett. 



A similar expansion holds for TrH n °^\ with coefficients a^o- As e — >• 0, Mazzeo and Rowlett 
show that a , e —> ao,o an d a i,e ~ ► a i,o, but that 



where the sum is taken over all interior angles oij of Qq |MR] . This anomaly in the heat trace 
asymptotics was first discovered by Fedosov [Fedj in a context unrelated to degeneration. Its 
expression was first simplified by Ray; Ray's work is unpublished, but a clear explanation 
may be found in a paper of van dan Berg and Srisatkunarajah [vBSJ. 

Mazzeo and Rowlett analyze the asymptotic structure of TrH ne (t) for small e and t and 
use this structure to explain the non-uniformity in the heat trace asymptotics. Their work 
motivates our approach to analyzing the determinant under conic degeneration; many of the 
features in this setting, including the non-uniformity in the asymptotics, are mirrored in our 
case. 

The behavior of the heat kernel on our degenerating family of manifolds has also been 
studied by Rowlett [Row], who outlines a direct construction for the heat kernel on Q t in a 
slightly more general setting known as 'asymptotically conic convergence.' This construction 
uses a variant of the analytic surgery technique of McDonald. Unfortunately, certain details 
in [Row] do not seem to be correct. One lesson from this construction is that the behavior 
of the full heat kernel under conic degeneration is quite complicated. In the present work, 
we introduce a less involved approach which focuses only on the heat trace. Since the zeta 
function only depends on the heat trace and not on the off-diagonal heat kernel, our approach 
suffices for the analysis of the determinant. 

1.3. Main results. Many of the proofs in the present work depend on the asymptotic 
structure of the heat kernel on Z, which is analyzed in |Slj . In particular, we prove in |S1] 
that it is possible to define a renormalized zeta function and determinant of the Laplacian 
on Z, denoted R (z and R det Az respectively. We then have the following approximation 
formula for the determinant under conic degeneration: 
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Figure 3. The blown- up space Qq. 

Theorem 2. As e ^ 0, 

log det Aq e = -2 log e( R ( z (0)) + log det A Qo + \og R det A z + o(l). 

To prove Theorem [21 we use the representation of the zeta function in terms of the heat 
trace, and then study the precise microlocal structure of the heat trace on Q t as a function 
of t and e. We do not use Rowlett's analytic surgery approach; rather, we perform a direct 
parametrix construction for the heat kernel of Q e , which gives us enough information to 
pass to the heat trace. The key structure theorem we prove is motivated by the work of 
Mazzeo and Rowlett in [MRj . Let Q be the quadrant x K+(e), and then let Qo be 

[Q; {(0,0)}]; that is, Q with a radial blow-up at (0,0) in the coordinates (y/t,e). The space 
Qo is illustrated in Figure EJ and its boundary faces are labeled L, F, and R as illustrated. 
Let Xi be a smooth radial cutoff function on Q e , equal to 1 on r < 15/16 and on r > 17/16, 
and let Xi = 1 ~~ Xi- Then we have the following: 

Theorem 3. TrH n "(i) is polyhomogeneous conormal on Qo, with leading orders 
and F (in terms of y/t) and leading order at R. Moreover, we have 

TrH Qe (t) — [ xi(^)H z (^,z,z) dz + / X 2{z)H Q °{t, z, z) dz + R(e,t), 

where R(e, t) is polyhomogeneous conormal on Q with infinite- order decay at t = 

Theorem [3] both describes the asymptotic structure of the heat trace and identifies all the 
terms in the asymptotics at L and F. On the other hand, it does not identify the leading- 
order term in the asymptotics of the heat trace at R, which is necessary for the proof of 
Theorem [2j We therefore need the following heat trace convergence result as well: 

Theorem 4. As e — > 0, for each fixed positive t, 

TrH n '(t) -> TrH n °(t). 

Moreover, as t — > oo, TrH Qf -(t) — 1 decays exponentially in t, uniformly in e. 

The natural follow-up to Theorem [2] would be to obtain a similar determinant approx- 
imation formula for the Laplacian acting on differential forms, preferably also on sections 
of twisted form bundles. This would enable us to derive an approximation formula for the 



— n at L 
(3) 



6 



DAVID A. SHER 



analytic torsion. If we had a similar approximation formula for the Reidemeister torsion, we 
could hope to apply the Cheeger-Muller theorem on each Q e and then take a limit as e — > 
to obtain a Cheeger-Muller theorem for Qq- This is ongoing work. 

1.4. Outline of the proof. The proofs proceed in two steps: first we assume Theorem [3] 
and prove the rest of the results, and then we return to prove Theorem [3j In section 2, we 
extend Theorem [TJ to cover dimensions 2, 3, and 4. We then prove a uniform lower bound 
on the eigenvalues of A^ e . Combining these two results with Theorem [3] allows us to prove 
Theorem HJ Then, in section 3, we use Theorem [3] and Theorem [4] to prove Theorem [2J We 
do this directly, by using the heat trace to construct the zeta function and then analyzing 
the derivative at zero. This analysis is an example of how geometric microlocal structure 
theorems such as Theorem [3] may be powerfully combined with geometric information such 
as Theorems [TJ and HI 

Section 4 is devoted to the proof of Theorem [3j We construct a parametrix for the heat 
equation on Q € by smoothly gluing together the heat kernels on eZ for r < 1 and on Qq 
for r > 1. The heat kernel on eZ is a parabolic scaling of the heat kernel on Z, which 
explains why we need to understand the long-time heat kernel on Z as well as the short- 
time heat kernel. We combine the results in |Slj describing the structure of the heat kernel 
on Z with the results in |Moj on the heat kernel on manifolds with conic singularities to 
prove Theorem [31 This section involves extensive use of geometric microlocal analysis, in 
particular Melrose's pushforward theorem, and is the most technical part of this work. For 
the necessary background on geometric microlocal analysis, see |Me| IMe2| IGr| (Sj [ST|. The 
Appendix contains the proof of one particularly technical result. 

1.5. Acknowledgements. This work comprises the second part of my Stanford Ph.D. thesis 
[S]. First and foremost, I would like to thank my advisor, Rafe Mazzeo, who introduced me 
to this problem and shared tremendous advice and support. Among the many other people 
who provided insight and helpful conversation, I would like to particularly mention Xianzhe 
Dai, Leonid Friedlander, Colin Guillarmou, Andrew Hassell, and Julie Rowlett. Finally, I 
am grateful to Gilles Carron and the Universite de Nantes for their hospitality during fall 
2010, and to the ARCS foundation for their financial support during the 2011-2012 academic 
year. 

2. Heat trace convergence 

In this section, we prove Theorem [U Our strategy is to combine the spectral convergence 
of |Ma2j . [Row], and |AT] with a uniform lower bound on the eigenvalues. The uniform lower 
bound will allow us to control the remainders in the infinite sum 

oo 

TrH a '(t) = ^e~ K ' kt . 

k=0 

2.1. Spectral convergence. As discussed in |Row2] . the spectral convergence result of 
Theorem [TJ has only been demonstrated when n > 5. We must therefore extend Theorem [TJ 
to the cases where 2 < n < 4. Our approach is the same as that in jMa2] and [R owl . The 
key is to prove that if we have a sequence of eigenvalues X e of Aq e which converge to A, then 
A is an eigenvalue of the Friedrichs extension of Aq . This will show that every accumulation 
point of cr(A ne ) is in cr(A no ). 
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Assume that for each e > 0, f e is a nonzero eigenfunction of A e with eigenvalue A e , 
normalized so that ||/ e ||oo — 1- Over any compact set in Q away from the conic tip, the 
Arzela-Ascoli theorem implies that f e converges uniformly to a limiting function /. By 
elliptic estimates, the convergence is actually in C°°, so the limit / is smooth and satisfies 
Aq / = A/. Note that since each fj is bounded by 1, we also have ||/||oo < 1> which 
automatically implies that / is in the Friedrichs domain of Aq . In order to show that A is 
an eigenvalue of Aq , we must show that / is nonzero. 

Following |Ma2] and |Row] , we define a weight function R e on Q e for each e. We set R t = e 
for r < e, R e = r for e < r < 1, and R € = 1 outside r = 1. Let 5 be a small positive number 
which we will choose later in the argument. For each e, let p e be the point at which the 
supremum of |/ € -Rf| is attained. Re-normalize the f e by writing ip e = /e|-Rf/ e (p e )| _1 , so that 
the supremum of |-Rfy e | is 1 and it is achieved at p e . Note that the ip e are still eigenfunctions 
ofA^ 

We now analyze the behavior of R € {Pe) as e approaches zero. Passing to a subsequence, 
we may assume that we are in one of three cases. 

Case 1: Suppose that R e (Pe) approaches R > as e goes to zero. In this case, some 
subsequence of p 6 approaches some p in Q away from the conic tip. Arguing as before, the 
(f e converge to some function (p on compact subsets away from the conic tip, and since the 
renormalizing factor R s e (p e ) converges to R s , we have (p = R~ s f. Since \R^<p e {p e )\ = 1 for 
all e, we must have \<p{p)\ = R' 1 - Therefore \f(p)\ — 1, so / is nonzero; we already showed 
that / is in the Friedrichs domain, so we conclude that A G ct(Aq ). 

Case 2: Suppose that there exists a constant C so that R € {p € ) < Ce as e goes to zero. 
We will derive a contradiction. Restricting ip e to the region where r < 1, we consider ip e as 
an eigenfunction of A eZ with eigenvalue A e . Rescaling by e, we view (p e as a function on Z; 
for \z\ < 1/e, it is an eigenfunction of Az with eigenvalue e 2 A e . 

For each e, let ij) e = (R t p t ) S (p e ; it is also an eigenfunction of A z with the same eigenvalue, 
and we have | -0 e (^ e ) | = 1. Moreover, p e maximizes R^\ip e \, so it also maximizes R & t \ip e \. Since 
R € rescales to a function on Z which is e\z\ outside \z\ = 1, we have that for \z\ > 1, 

(e\z\) s \A(z)\<Re(Pe) S <C s e 5 . 

Therefore, for \z\ > 1, we have (V'el- 2 )! ^ C" 5 ^! -5 . 

As before, we take a limit of the ip e (uniformly on compact subsets of Z) to obtain a limit 
ip. Since e 2 A e goes to zero, ip is a harmonic function on Z. The rescaled points p e lie in a 
bounded subset of Z and hence converge to some p G Z, and so we must have ip(p) = 1. 
Since each |V'e( 2; )| — C s \z\~ s for \z\ > 1, the limit ip must satisfy the same bound. However, 
since 5 > 0, ijj is a nonzero harmonic function on Z which approaches zero at infinity. This 
is impossible by the maximum principle, and we have a contradiction. 

Case 3: Finally, suppose that R e (p e ) approaches zero but that e _1 i? e (p e ) goes to infinity. 
Consider the restriction of ip e to the region in Q e between e and 1, where the metric is exactly 
conic and R e is precisely r. We again rescale, but this time by R e (pe) instead of e. The result 
is that we view (p t as a function on the region of the infinite cone Cn between s = e/(R e (p e )) 
and s = l/(R e (p e )), where s is the radial variable on the cone. Note that p e has s = 1. 

As in Case 2, we again let ip e = R e (p e ) S (p t . This rescaled function ip e is an eigenfunction 
of Ac N with eigenvalue R e (p e ) 2 X e , and \if> e \ attains its maximum of 1 at p e . The function 
R e lifts to R e (p e )s, so we have that for any point (s, 0) G CV with s between e/(R e (p e )) and 
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(R € (p € )s) s Ms,6)\<R € (p € ) s . 

Hence |^ e (s,0)| < s~ s . 

As e goes to zero, the lower bound for s goes to zero and the upper bound for s goes to 
infinity. Therefore, as in the other cases, the functions ip e converge to a limit ip on C^. The 
rescaled points p € all have s — 1 and hence have a subsequence converging to a point p. 
The eigenvalues R € (p € ) 2 X € converge to zero, so ip is a harmonic function on Cn', moreover, 
■0(p) = 1, so it is nonzero. Finally, ip(s,9) < s~ s for all s. If 6 is not an indicial root of 
the conic Laplacian, this is impossible (see [Ma2j). Since the indicial roots are discrete, we 
may pick 5 which is not an indicial root. This contradiction completes the proof that A is in 

In order to complete the proof of spectral convergence, we must also show that every 
point of er(An ) is an accumulation point of cr(AnJ, and that the multiplicities are correct. 
However, the proof of these assertions in [Row] does not depend on the assumption n > 5. 
Therefore, the remainder of the proof proceeds as in [Row] . and we will not repeat it here. 
We have now established Theorem [T] for all n > 2. 

2.2. Lower bounds for eigenvalues. We now prove uniform lower bounds on the eigen- 
values of Aq e . This sort of lower bound is generally difficult to prove using classical methods. 
However, as we will see, it is a relatively straightforward consequence of Theorem [3j The 
key step is the following upper bound on the heat trace: 

Lemma 5. For any fixed T > 0, there is a constant C so that for all t < T and all e < 1/2 
(including e = 0), 

TrH Qe (t) < Cr n/2 . 

Proof. From Theorem |3l TrH ne (t) is polyhomogeneous conormal on Qo, with leading orders 
— n at L and F (in terms of yi) and at R. The function t n l 2 is polyhomogeneous on Q and 
hence on Qq, with leading orders n at L and F and at R. Therefore, t n ^ 2 TrH nE (t) is poly- 
homogeneous on Qo with leading order at each boundary face. In particular, t n ^ 2 TrH nE (t) 
is continuous up to the boundary of Qo. If we restrict to t < T for any fixed T and to 
e < 1/2, then t n / 2 (Trif° e (t)) is uniformly bounded, proving the lemma. □ 

Lemma immediately implies a uniform lower bound for the eigenvalues: 

Lemma 6. There is an Nq G N and a C > such that for all k > Nq and all e < 1/2 
(including e = 0), 

Xe,k > C'0 n . 

Proof. Let N e (X) be the number of eigenvalues of A^ E less than or equal to A. Then apply 
Lemma [5] with T = n/2 to obtain that for all t < n/2, 

N e (X)e- xt < e ' Kjt ^ TrH Q <(t) < Cr n/2 . (4) 

A £ ,i<A 

Therefore, for all t < n/2, 

N e (X) < Ct- n / 2 e xt . (5) 

The right side of (jSJ) is minimized when t = (n/2)/X, which is bounded by n/2 when A > 1. 
So when A > 1, 

N t (X) < CX n/2 . (6) 
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Applying this to A = A e> fc, we see that for A £) fc > 1, 

k < ck[1 

and hence 

K, k > c'k 2 ' n . 

Finally, note that N e (l) < C, so picking any N > C ensures that if k > N, A e ,fc > 1. This 
completes the proof of Lemma [6j □ 

This section illustrates the power of microlocal structure theorems such as Theorem [31 it 
is possible to prove Lemma [6] via a classical argument involving isoperimetric inequalities, 
but that method is much more involved. 

2.3. Heat trace convergence. We now prove Theorem HI starting with the exponential 
decay claim. By Theorem!]], A £) i converges to Ao,i > 0, and hence is bounded below uniformly 
by Ao,i/2 for small e; the same is of course true for for all k. Combining this observation 
with Lemma O we have 

oo 

TrH^{t) < 1 + (N - l) e -( A o, 1 /2) i + ^ e -*7'**/*_ (7) 

fc=JV 

Since the infinite sum in ([7]) is convergent for any t and exponentially decaying in t, sub- 
tracting 1 from both sides immediately proves the second statement in Theorem HI 

To prove the heat trace convergence part of Theorem HI fix 5 > and t > 0. For any 
M > Nq, we write 

M oo oo 

\TrH Qc (t) - TrH n °{t)\ < ^ |e" A - fe * - e" AE -°'| + ^ e^ 1 + e" A °' fc *. 

k=l k=M+l k=M+l 

Lemma [6] implies that the second and third terms are both bounded by 



-tC'k 2 / n 

& 1 



E 

k=M+l 

which goes to zero as M goes to infinity. So we may fix an M such that both of these terms 
are bounded by 5/3 for any e < 1/2. Finally, by spectral convergence, for sufficiently small 
e, the first term is also bounded by 5/3. We conclude that for sufficiently small e, 

\TrH nc (t) -TrH n °(t)\ < 5. 

Since d > was arbitrary, this completes the proof of Theorem HI 

3. The determinant formula 

The proof of Theorem [2] involves direct analysis of the zeta function on Q e , in particular 
the standard heat trace formula 

CnAs) = W ~, / (TrH e (t) — l)t s-1 ds, (8) 

r(s) Jo 

where equality is in the sense of meromorphic continuation from the region where 3?s > n/2. 
We prove Theorem [2] by decomposing the heat trace into several pieces and analyzing (JSJ) 
for each piece, using Theorems [3] and @] 
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3.1. Long-time contribution. We begin by studying the long-time piece of the zeta func- 
tion: 



1 f°° 



ds. (9) 



The integrand decays exponentially in t as t — > oo for any s G C, so we can differentiate 
under the integral sign. The function — Sjy is holomorphic at s = with leading order s, so 
differentiation gives: 



/oo 
(TrH^(t) - I)*" 1 dt. 



(10) 



As an aside, note that in fact (q e l (0) = 0; the integral in ([9]) is holomorphic on C, and 
hence multiplication by the inverse gamma function gives a holomorphic function which is 
at s = 0. 

By Theorem HI there exist /i > and C < oo independent of e such that: 

\{TrH n ^{t) -1)| < Ce""*. 

Moreover, by the same theorem, 

limTrH^it) - 1 = TrH n °(t) - 1. 

We may therefore apply the dominated convergence theorem to (TTUT) to conclude that 

hmC^, L (0) = ^ o , L (0). (11) 

In fact, we expect that there should be a polyhomogeneous expansion for L (0) as e goes 
to zero, but it is difficult to analyze the long-time behavior of e-derivatives of TrH e (t). 

3.2. Projection term. Now we consider the short-time zeta function 

Cn e ,s(s) = J-j JjTrH^t) - l)t s_1 cfo. (12) 
By direct computation, we may write 

( ^ {8) = ikl! TrH "' {t)f " d8 -^y (13) 

The second term in fflBl is the contribution from the projection off the constants. 

3.3. Leading order term at R. Now that we have analyzed the long-time and projection 
contributions, we consider the integral 



1 r 1 

— / TrH a '(t) t 8 ' 1 dt. 
(«) Jo 



(14) 



. To best take advantage of our knowledge about the leading-order terms of TrH ne (t) on Q , 
we break off the leading-order term at R in a neighborhood of R. In particular, fix b G (0, oo), 
and break the analysis of ( fl4l) into two regions: \/t < be and y/i > be. We work in local 
coordinates in each region; let r = t/e 2 and r\ = e/y/i, as inspired by [MR] and illustrated 
in Figure [3J When \/t < be, we use the coordinates (y/r,e); when \/t > be, we instead use 
(v,Vi). 
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In the region {\/t > be}, we break off the leading-order term at R, which by Theorem H] is 
TrH n °(t)x{Vi > be}, and consider it separately. Its contribution to (TH|) is precisely 

-i- T TrH^ity- 1 dt = (Cn ,s(s) + -^) - =^r f ' TrH^{t)t s - 1 dt. (15) 

Note that (|15j) is initially valid for large s; however, since both sides have meromorphic 
continuations, it is valid for all s. 

As proven by Cheeger [Ch2] and generalized by Bruning and Seeley |BSJ, TrH Q °(t) has a 
short-time heat expansion given by 

n 

TrH n " (t) = a k t {k - n)/2 + K log t + TZ(t) , 

k=0 

where TZ(t) is 0(t a ) for some a > 0. We plug this expansion into fll5p . All the integrals 
except for the error term may be evaluated directly, and yield: 

-fn E^Tr^wrr- - ^frr (i - (21og6 + 21oge)) - 

r ( s )^ (/c-n)/2 + s sr(s) s 

This is a meromorphic function on C, with a simple pole at if K ^ 0. 

As for the error term, it is of order t a for a > 0, so its contribution to f JT5|) converges and 
is holomorphic in a neighborhood of s = 0. It is bounded by 

Ct a+S ^ dt = —— ^ -(b 2 e 2 ) a+s . 



T(s)J T(s)(a + s) 

In the neighborhood \s\ < a/2, this is a holomorphic function for each e which decays in 
e. Thus the contribution of this error term to the zeta function is o(l) in e. To show that 
the same is true for the determinant, we take an s-derivative of the R(t) portion of (TT5|) and 
obtain: 

pb 2 e 2 i rb 2 e 2 

- Rity- 1 dt - —- RitU'-Hogt dt. (16) 
Jo T ( s ) Jo 

Essentially the same analysis, using the fact that R(t) < Ct a , shows that (TT6|) is o(l) in e. 

Combining the contributions of the long-time zeta function, the projection off the con- 
stants, and the leading-order term in \/i > be, we have: 

1 n (U F \k-n+2s Kb 2s f 2s 1 

<*W - W) ^ ( L' n)l2 + s ~ ir^T'j - ^ + 2i ose)) , (17) 

plus a term whose contribution to the determinant is o(l) as e goes to zero. 

3.4. The Q Q term. Now we consider the portion of ( Tl4|) which remains after subtracting 
the leading-order term at R. Write 

TrH n <(t)= [ X i{z)H ez {t,z,z) dz + [ X 2{z)H Q °{t, z, z) dz + R{e,t). (18) 

We analyze these three terms' contributions to ( I14p separately; however, we must subtract 
the leading order term at R in the region y/t > be. By Theorem [3] and its proof, each of the 
three terms is phg conormal on Qq with leading order at R, so we may separately subtract 
the leading order term at R from each term in (TTJ 
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First, we consider the f2 term. Its contribution to (!T4l) is 

1 



r( s ) 



(/ X 2(z)H nQ {t,z,z) dz)^ 1 dt; 
o Jn 



however, we have to subtract the leading-order terms at R in y/t > be. The Qq term is 
independent of e, so the leading-order term at R is the whole integral. Therefore, we only 
integrate from to b 2 e 2 . The contribution to (TL4]) is 

b 2 e 2 

X2 (z)H Qo (t,z,z) dz)t s ~ l dt. 

n 



T(s) 

This analysis is similar to the previous section. Note that the short-time heat expansion for 
f^o is local, and X2{z) localizes away from the cone point, so there is no logarithmic term. 
Therefore, there exist dk such that 



„ n 

j X 2(z)H n °(t, z, z) dz = J2 d k t (k - n)/2 + C(t 



As in the discussion of the leading-order term at R, we obtain a contribution to (j!4j) of 

1 - (hAk-n+2s 

fnE S* ( » (19) 

r ( s )^ (k-n)/2 + s 
modulo an error term whose contribution to the determinant is o(l) in e. 



3.5. The error term. Now we consider the term coming from R(e, t). In the region y/t < be, 
this decays to infinite order in both r and e; in particular, it is less than Cre 2 = Ct for some 
fixed C. As before, any term of order 0(t a ) for a > contributes only o(l) to the determinant 
when integrated over the region \ft < be, so the contribution from this region is o(l). 

On the other hand, where \/i > be, R(e, t) is a function of (77, y/t). It decays to infinite order 
in y/t, and since we have subtracted the leading order term at R it decays to some positive 

3(5 

order 5 in 77. In particular, R(e, t) is bounded in the region xft > be by Cy/t 7] s = Ct 5 e 5 . 
The contribution to fTl4l) is an integral between t = b 2 e 2 and t = 1, so it is convergent and 
holomorphic for all s (we do not integrate down to s = 0). We may then differentiate under 
the integral sign and obtain that its contribution to 5(0) is 

f R(e,t)r l dt, 

which is bounded by 

Ce 5 Ct^dt = %e s . 
Jo 

Therefore, the remainder term R(e,t) contributes only o(l) to the determinant as e ap- 
proaches zero after we subtract the leading order at R. 
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3.6. The eZ term. It remains only to analyze the contribution to (jl4]) of the first term in 
(TI8|) . When we plug this term into ( TH| . we get 

(/ X i{z)H eZ {t,z,z) dz)^' 1 dt, 



which may be rewritten, by scaling and a change of variables, as: 

2e 2s rl/e 



We now recognize the term 



f \ f Xi(ez)H z (r,z,z)dz) v ^ 2S ' 1 dyfF. 
Jo J z 



Xi(ez)H z (r,z,z) dz (20) 



from the discussion in |Slj of the renormalized heat trace on asymptotically conic mani- 
folds. In [51] . we prove that (1201) is polyhomogeneous conormal in (r, e) for small r and 
on X^r- 1 / 2 ^) for large rfl Moreover, since Z is exactly conic near infinity, [SI] gives an 
explicit expansion for ( 1201) as e goes to zero: 

„ 71— 1 

/ XiMH Z {r,z,z) dz = Y,D k r^l 2 t k - n - Cloge 

J Z k=0 

+ ( R TrH z (r) - Cl log ) + R(e, r), (21) 

where R(e,r) vanishes as e — > for each r. Moreover, if we let u n (z) be the coefficients of 
the heat asymptotics on the infinite cone CV, then 

C = -[ u n (l,y) dy, D k = — — / u k (l,y) dy, 
J n K — n J N 

2 /-2 



l k = -l x'i( r ) r n dr,l log = - ^(r) logr dr. 

Jl/2 J 1/2 

Note that by the phg conormality of (T2Q|) and the expansion (1211) . the remainder term R(e, r) 
has the same phg conormality as (TSUI) . 

We will plug (j2"T]) into Since we still need to subtract off the leading order terms in 
y/r > b, we break up the integral at y/r = b. 

3.6.1. t < b 2 : For r < b 2 , (e, r) are good coordinates, and we rewrite the integral ( 114]) as an 
integral in y/r, from -y/r = to y/r = b: 

2e 2s 



Plugging the expansion ( |2T]) into ( |22l) gives 

n-l 

E 

fc=0 



/ TrH^{T 2 e)s/T S ^d^. 
Jo 



(22) 



2D k e 2s+k - n b 2s+k ~ n Ce 2s b 2s log e Cl log e 2s b 2s 
(2s + k-n)T(s) sr(s) sf(s) 



1 The definition of A^fY" 1 / 2 , e) may be found in [MS]; it is just M + (t _1/2 ) x M+(e) with a radial blow-up 
at t -1 / 2 = e = 0. As always, see |Me| IMe2) IGr] for background on blow-ups and polyhomogeneous conormal 
functions. 
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f fA f R TrH z {r)^ 23 - 1 dyfF + R(e,s) : 
r 0) Jo 



>o 

where 

R(e, a) = ^r f R(e, r)^ 1 dyfP. (23) 
1 V s ) Jo 

Note that R(e, r) is phg conormal in (e, r) for small r and the leading order at e = is 
positive, so the coefficients of powers of r all decay to a positive order in e; let 5 be less than 
that positive power. Then there exist I, z^, and for each k between and I such that 



R(e, T ) = ^6 fc (e)r^(logr) Pfc + C(eV/ 2, 

k=0 



where fefe(e) are 0(e 6 ) as well. Plugging this expansion into (I2"3"]h the 0(e 5 r 5 ^ 2 ) term is 
bounded by CeV/ 2 = t s/2 for y/r < b, and hence its contribution to the zeta function and 
determinant is o(l), as before. The rest of the expansion gives 

7^X>( e ) f r z ^ 2s - l {\ogr)^ d^. 

Each term, in a neighborhood of s = 0, is e 2s+s times an explicit meromorphic function of 
s with coefficient depending on b. Therefore, each of the coefficients in the Laurent series 
for (1231) around s = is 0(e 5 ), and so the contribution of this remainder term to the zeta 
function and determinant is o(l) (in fact, 0(e s )). 

3.6.2. r > b 2 : In this region, we use the coordinates (77, and the expansion ( |2T]) becomes: 

n— 1 

£ ^(V^*-" - C log( V / t) - C log 77 + K TrH z ( V - 2 ) - Cl log + R. (24) 

fc=0 

We need to subtract the leading order term at R, which in these coordinates is the face 
77 = 0. We do this first for the non-remainder terms and then for R. 

The expansion ([21]) is phg conormal on Q , hence in (77, \/t) in this region, and it has 
leading order at R. The non-remainder terms are phg conormal as a consequence of the 
phg conormality of the renormalized heat trace (proven in [SI]). Moreover, they are precisely 
the finite and divergent parts of (|2Q|) at F, and hence have leading order at R; this implies 
that R TrH z (t]~ 2 ) has leading-order behavior of Clog 77 at 77 = 0. If we let be the constant 
term in the expansion at rj = of R TrH z (r]~ 2 ), we see that the terms which remain after 
subtracting the leading order at R of the non-remainder terms in (|24|) are precisely 

- Clog 77 + R TrH z ( V - 2 ) - /«,. (25) 

We now plug these terms into (JHJ); we could switch (JH|) to an integral in 77, but it is easier 
to switch (1251) back to (e,r), and it becomes 



R t 1 ^ uZ I 



Clog(v^) + R TrH z (r)-f OQ . 
Plugging these into (122]) and integrating from y/r = b to yfr — 1/e gives a contribution of 

2e 2s f l/e 



J Vlog(v^) + R TrH z (r) - U^'' dy/t. (26) 
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To get rid of the artificial upper bound of 1/e, examine 
2e 2s 



T(s) 



poo 

/ (Clog(v^) + R TrH z (r) - f^y/^ 1 dy/¥. (27) 

Jl/e 

Because of the positive-order decay of ( 1251) at r = oo, the integral in ( 1271) is well-defined in a 
neighborhood of s = and has uniform-in-s decay to a positive order in e. The same is true 
of its s-derivative, precisely as in subsection 13.31 Therefore, ( 1271) contributes only terms of 
order o(l) to the determinant, and we can replace the analysis of (|26|) by analysis of 

2e 2s 
T(s) 

We can directly evaluate most of the integrals in (128 p . and we obtain a contribution of 



POO 

/ (Clog(v^) + R TrH z (r) - f^y/^ 1 dy/¥. (28) 
Jb 



Ce 2s 6 2 * , , 1. 2e 2s f°° - ttZ/ n ^s-i , r- fooe 2s b 2s 



2sT(s)^ b '^ + W)l RT ^(r)^-dV-r + ^. (29) 

As for the remainder, since both ( I24p and the non-remainder terms in (|240 are phg conor- 
mal on Q , so is R; moreover, it has positive leading order at F and leading order at R. 
When we subtract the constant term at R, R has positive-order decay at both L and R, 
hence at both rj = and \/t = 0. As in the analysis of the remainder term in subsection 13.5} 
R only contributes a term of o(l) to the zeta function and determinant. 

Combining all the analysis in this subsection, we see that the total contribution of the eZ 
term to the zeta function, after subtracting the leading orders at R, is 

g 2D k e 2s+k - n b 2s+k ~ n Ce 2s b 2s log e Cl log e 2s b 2s 



k=0 



(2s + k-n)T(s) sT(s) sT(s) 



-[2logb--] + ^-^- + 4-- I "TrH^r)^- 1 dr, (30) 



2sr(s) s s sT(s) T(s) J 

modulo terms which contribute only o(l) to the determinant expansion. We now recognize 
the last term above as precisely € 2s ( R (z(s)) - that is, e 2s times the renormalized zeta function 
on Z. 

3.7. Proof of Theorem [2J, We now combine (fTTj) . (jT9j) . and (1301) . treating the k = n case 
of (|T7|) . ( !T9|) . and (|30|) separately. The zeta function Co e (s), modulo terms which contribute 
only o(l) to the determinant, is: 

1 yi (-a* + a k + D k )e 2s+k - n b 2s+k - n (C + 2K)e 2s b 2s 
^) + —^ a + (fc _ n) /2 + 2^r( s ) 

+ ^(-a n + ~a n - C7 log - (C + 2K)(logb + loge) + /«,) + e 2s (%(s)). 

This expression is a meromorphic function of e. We may therefore evaluate the s coefficient 
of its Laurent series at s = 0. When we do this, we must take derivatives of the gamma 
function, so we write the Laurent series for as s + A^s 2 + A 3 s 3 + . . .. We obtain 

C^ E (0) = CnM + R Cz(0)(2 loge) + R <* (0) + £ ^ + g + 

£^ (k-n)/2 
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+ (4» + A 2 (2 log b + 2 log e) + ~ (2 log b + 2 log e) 2 ) 

+(2 log b + 2 log e + A 2 ) (-a n + S n - C7 log - (C + 2K) (log b + log e) + U) + o(l) . 

This equation appears quite complicated. However, the choice of 6 was arbitrary, so the 
coefficients of all powers of e must be independent of b. This allows successive simplifications. 

First notice that the coefficient of e k ~ n , for each k between and n — 1, is a constant times 
b k ~ n . This coefficient must be independent of b, so the constant must be zero. Therefore 
— afc + <2fc + .Dfc = for all k G [0, n — 1]. Next, we isolate the coefficient of (loge) 2 ; it is 

° + 2K ■ 2 - 2(C + 2if) = -(C + 2K), 

which is independent of b. However, the coefficient of log e is 

-2(C + 2*0 logfe + 2(-a n + a n - Cl log + /«, + R Cz(0)). 

This coefficient must be independent of 6, which forces C + 2K = 0. Using the two identities 
we have discovered, we simplify the expression for Cn £ (0) an d obtain 

Cn o (0) + R Cz(0)(21oge) + R Cz(0) + (21og6 + 21oge + A 2 )(-a n + a n - Cl log + /„), 

plus terms vanishing as e goes to zeroB Finally, the constant term in this expansion is 

Cn o (0) + R Cz(0) + (21og6 + A 2 )(-a„ + a n - C/ log + /«,), 

and again it must be independent of 6. Therefore — a n + a n — C/i og + foo = 0. We can now 
write 

C^ £ (0) = 21oge(%(0)) + C^ o (0) + R Cz(0) + o(l). 

Multiplying by —1 and using the definition of the determinant of the Laplacian completes 
the proof of Theorem [21 

As a quick check on Theorem [21 suppose that Z = R n ; this is the trivial case, where Qq 
is flat in a neighborhood of P, with no conic singularity, and Q e = Qq for every e. Then it is 
easy to compute the renormalized heat trace; it is the finite part as e — » of 

/ H^(t,z,z) dz = V ° l ( ( "^! e) = e-^t)-^Vol(r < 1). 

Jr<l/e {QITZ) I 

The finite part is precisely zero, which means that the renormalized heat trace is identically 
zero, and hence the same is true for the renormalized zeta function and determinant. The- 
orem [2] then reduces to Cn E (0) = Cc o (0) + °(1)> as expected (in fact, the error is zero in this 
case). 



4. The parametrix construction 

In this section, we construct the heat kernel for Q e by gluing together the heat kernels on 
eZ and Q . We then use the construction to analyze the fine structure of the heat trace and 
prove Theorem [31 



These cancellations may seem rather fortunate, but they just reflect the artificial nature of the decom- 
position at y/r = b. 



CONIC DEGENERATION AND THE DETERMINANT OF THE LAPLACIAN 



17 




Figure 4. Cutoffs as functions of r. 



4.1. Outline of the parametrix construction. In order to do a gluing construction, we 
need to introduce cutoff functions. First restrict attention to < e < 1/2, and note that 
Q e is exactly conic between r = | and r = 2. Recall from the introduction that xi is a 
smooth radial cutoff function on Q t , equal to 1 on r < 15/16 and on r > 17/16, and that 
X2 = 1 — Xi- Define two slightly larger radial cutoff functions; let xi be 1 on r < 9/8 and 
on r > 5/4. Let X2 be on r < 3/4 and 1 on r > 5/8. These cutoffs are illustrated, as 
functions of r, in Figure 6.1. 

We can now define the initial parametrix: 

Definition. The parametrix G e (t) is the t- dependent operator on Q e with Schwarz kernel 

G e (t, z, z') = Hz)H eZ (t, z, z')xi(z') + X2(z)H n °(t, z, z') X 2(z'). (31) 

Our construction is similar to the usual construction of the heat kernel, which is explained, 
for example, in [Rj. Denote by * the operation of t-convolution: if two time-dependent oper- 
ators A and B on a manifold M have Schwarz kernels A(t, z, z') and B(t, z, z') respectively, 
and the integral 



/ / A{t — s, z, z")B(s, z", z') dz" ds. 
Jo J M 



(32) 



converges for all t, then the t-convolution A * B is well-defined with Schwarz kernel given 
by (!32l) . It is easy to check that t-convolution is associative; denote by A* k the operator A, 
i-convolved with itself k times. It is also easy to check that Tr(A * B)(t) = Tr(B * A)(t) 
whenever both A* B and B * A are defined. 

We first define the error E £ (t, Z j Z ^ clS cL kernel on fi f : 

E e (t, z, z 1 ) = (d t + A Q( , z )G e (t, z, z'). (33) 

Then define 

oo 

K e (t, Z ,z') = Y,(-l) k+1 E: k (t,z,z'). (34) 

k=l 

We claim: 

Proposition 7. For each e, H Qe (t, z, z') = G e (t, z, z') - (G e * K e )(t, z, z'). 

Proof. The proof of Proposition [7] is relatively standard, and is based on [R]. We first notice 
that for any smooth function f(z') on Q e , 



hm / G e (t, z, z')f(z') dz' = f(z). 
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This is an immediate consequence of ( )3T|) and the definition of the heat kernels on eZ and Qo. 
Further, we claim that for each fixed e and each n, there is a C n e so that \A e (t, z, z')\ < C n ^t n 
for small t, where A is any of E, K, or G * K; this is a consequence of Lemmas [11] and 
O Assuming these two claims, and that all convolutions are well-defined, we compute 
(suppressing e in the notation): 

(d t + A Z )(G - G * K)(t, z, z') = E(t, z, z') - (d t + A Z )(G * K)(t, z, z') 
= E(t, z, z') -Urn [ G(t- s, z, z")K(s, z" , z') dz' - (E * K) (t, z, z') 

= E(t, z, z') - K{t, z, z') + (E* K)(t, z, z') = 0. (35) 
Moreover, for any smooth function f(z') on Q e , 

lim / (G - G * K)(t, z, z')f(z') = f(z) - lim(G * K)(t, z, z')f(z') dz' = f(z). (36) 

Uniqueness of the heat kernel completes the proof of Proposition [TJ modulo the proofs of 
Lemmas [11] and [HI □ 

Proposition [7] allows us to write 

TrH e (t) = TrG e (t) - Tr(G e * K e ){t) = TrG e (t) - Tr(K e * G e )(t). 

In fact, we notice that K e = E € — (K e * E e ), and so 

TrH e (t) = Tr{G € (t)) - Tr((E e * G e )(t)) + Tr((K e * (E e * G e ))(t)). (37) 

We will analyze each of the three terms and prove that each is polyhomogeneous conormal 
on Q ; ([3]) follows immediately. Throughout, we use the local coordinates on Q from the 
previous sections; see Figure [3] 

4.2. Heat kernels on eZ and Qq. We first need to understand the asymptotic structure 
of the heat kernels of eZ and f2 . Since f2 is a manifold with a single conic singularity, 
its heat kernel has been extensively studied. As discussed in section 3, there are short- 
time heat trace asymptotics due to Cheeger and Bruning-Seeley |Ch2t IBS] . However, for 
our parametrix construction, we need a complete description of the asymptotic structure of 
z,z') for short time; this description is provided by the work of Mooers [Mo] . Near 
the conic singularity, we use the polar coordinates z = (r, y) and z' = (r', y'), so that the heat 
kernel is a function on M 3 i _(\/t,r,r') x N y x N' y . Let M(t,r,r' ,y,y') be an iterated blow-up 
of this space: first blow up {r = r' = \/t = 0}, and then blow up the closure of the lift of 
= 0,r = r' > 0} (this process is illustrated in Figure [5]). Mooers then proved that: 

Proposition 8. |Mo] // Qq is any manifold with a conic singularity, and (r, y) are polar 
coordinates on Qq near the singularity, then H n °(t,r,r' ,y,y') is polyhomogeneous conormal 
on M(\/t,r,r' ,y,y') . The leading orders are —n (in terms of \/t) at the two blown-up faces, 
oo at the original \/i = face, and at the other two faces. 

As for the heat kernel on eZ, we need to understand H eZ (t, z, z') for bounded t. However, 
by the parabolic scaling of the heat kernel, we have 

H^(t,z,z') = e' n H z (- 2 ,-,-). 
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Figure 5. Asymptotic structure of the heat kernel on a conic manifold. 



Therefore, to understand H eZ for bounded time, uniformly as e goes to zero, we need to 
understand H z for both short and long time. The short-time structure is relatively well- 
understood as a result of the work of Albin [Alb] . and in |S1] we investigate the long-time 
structure by studying the low-energy resolvent on Z. 

We now summarize these results. On Z, let x = r -1 , so that a; is a 'boundary defining 
function' for infinity. Then the heat kernel is a function on (0, oo) t x x') x N y x N y >. 

Assume first that t < T for some fixed positive T; we construct the appropriate space in 
a series of steps. We form the 'scattering double space' Z^ c (x , y , x' , y') by starting with 
M? + (x,x') x N y x N y i and performing two blow-ups. First blow up {x = x' = 0}. Then 
let D be the 'lifted diagonal' in this blown-up space given by the closure of the lift of the 
interior diagonal {x = x' > 0, y = y'}, and blow up the intersection of D with the boundary 
x = x' = 0; the result is Z^ c (x,y,x',y'). Finally, form Z] c (x,y,x',y') x [0,T) t , and blow up 
the t = diagonal given by the closure of {t = 0, x = x' > 0, y = y'}. This space was first 
identified by Albin in |Alb] as the likely candidate for the polyhomogeneous structure of the 
short-time heat kernel on Z. In |Sf j . we prove as an immediate consequence of |Alb] that: 

Theorem 9. Fort < 1, the heat kernel H z (t,x,y,x',y') is polyhomogeneous conormal on 
[Z 2 sc (x, y, x', y') x [0, 1]^; {Vi = 0, x = x', y = y'}}. 

Moreover, there is infinite- order decay at all faces except the scattering face sc, where the 
leading order is 0, and the front face ff of the final blowup, where the leading order is y/i " ' . 
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Figure 6. The space Z\ % 



On the other hand, assume that t > T for some fixed positive T, and let w = t _1//2 , so 
that the heat kernel is a function of (w, x, x', y, y'). We define a blown-up space Z 2 SC as 
follows. First, begin with the space X^(w,x,x f ) x N' x N y r] is a 'b-stretched product' 
in the terminology of Melrose and Singer |MSj . More details may be found in |GHlj or 
[ST] . As before, let D be the closure of the lift of {x = x' > 0, y = y'}, then blow up the 
intersection of D with the face {w > 0, x = x' = 0} to obtain Z 2 sc . Z 2 sc has eight boundary 
hypersurfaces and is illustrated in Figure [6j the faces are also labeled as indicated in the 
figure. Let F(w,x,y,x',y') = H z (w~ 2 ,x,y,x',y'). In [SlJ, we prove the following. 

Theorem 10. For any n > 2, and for any fixed time T > 0, F(w, x, y, x', y') is polyhomo- 
geneous conormal on Z^ sc for t > T, where w = t~ l l 2 . The leading orders at the boundary 
hypersurfaces are at least at sc and n at each of bfo, rbo, Ibo, and zf, with infinite- order 
decay at lb, rb, and bf. 

Combining the short- and long-time structure theorems, we obtain a complete description 
of the asymptotic structure of the heat kernel on Z, illustrated in Figured we have indicated 
the leading order at each face (in terms of t 1//2 or w = t^ 1 ^ 2 ). In fact, as explained in [ST], 
both of these theorems hold even if Z is only asymptotically conic rather than exactly conic 
near infinity. 

Finally, we would like to be able to generalize the construction in this section to the 
case of the heat kernel for the Laplacian acting on sections of certain vector bundles over 
Z. In particular, we are interested in the twisted /c-form bundles needed to define the 
analytic torsion. Based on work of Guillarmou and Hassell in [GH2] and unpublished work 
of Guillarmou in [Guij . we expect that in many cases, the heat kernel for forms will have the 
same polyhomogeneity properties as above, but with different leading orders. In particular, 
the leading orders of the heat kernel for long time may be worse at rbo, Ibo, and zf. Therefore, 
we do not use the order n decay at those faces that we see in the case of the scalar Laplacian; 
instead, we use only the following strictly weaker assumption, which reflects the expected 
features of the problem for the twisted form bundles: 
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Figure 7. Asymptotic structure of the heat kernel on Z. 



Alternative Hypothesis. The heat kernel on Z has leading order at least at zf and at 
least n — a at rbo and lb , where a 6 [0, n/2). However, the coefficient of the term of order 
at zf decays to order strictly greater than n at bfo. 



4.3. The error terms E e and K e . To begin proving Proposition [7] (using only the alternate 
hypothesis), we first compute E e . Since H eZ and H n ° solve the heat equation on the regions 
r < 3/2 and r > 1/2 respectively, we obtain: 

E e (t,z,z') = iA^,xi(z)W eZ (t,z,z') X i(z') + [A Ue , z , x 2 (z)}H u °(t,z,z') X 2(z'). 

The commutators are first-order differential operators, with supports contained in the r- 
intervals [9/8,5/4] and [3/4,7/8] respectively. In these regions, Q e is exactly conic for all 
e < 1/2, so we may call the commutators V\(z) and V 2 {z). Then 

E e (t, z, z') = V 1 (z)H" z (t, z, z')xx(z') + V 2 (z)H n °(t, z, z') X 2{z'). (38) 

Note that V±(z) is supported in (9/8,5/4), while Xi( z ') is supported in r < 17/16; similarly, 
V 2 (z) is supported in (3/4,7/8), and X2( z ') is supported in r > 15/16. Therefore, the heat 
kernels in the definition of E e are each away from the diagonal. Off-diagonal heat kernels 
usually decay to infinite order at t = 0, which motivates the following lemma. Recall from 
section 2 that R t {z) is the weight function on Q e which is e inside r = e, 1 outside r = 1, 
and r in between; define a modification W e (z) to be r between r = e and r = 3/2 and 3/2 
outside r = 3/2. Then: 

Lemma 11. For each N > 0, there is a constant such that for all e < 1/2, 

sup \(W e (z')) a E e (t,z,z')\ <C N t N . 
z,z'en e 

Moreover, for each fixed z and z' both outside r = 1/2, E e (t,z,z') is phg conormal on Qo, 
and the dependence on z and z' is smooth. 

We require r > 1/2 in the second part of the statement in order to regard z and z' as 
coordinates on the fixed space Q . We now claim: 
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Corollary 12. For each fixed z and z' both outside r = 1/2, E e (t,z,z') is phg conormal on 
Q with infinite-order decay at \ft = 0. 

The corollary follows from Lemma [TT] and the general blow-up fact, from [Me2] . that any 
function phg conormal on [M; P] which vanishes to infinite order at the blown-up face is in 
fact phg conormal on M without the blow-up (see Proposition A. 3 of [5]). We now prove 
Lemma [TT1 

Proof. First consider the second term in E e , V2(z)H n °(t, z, z')x2( z ')- It is independent of e, 
and its support in (z, z') is away from the conic tip in both variables and also away from 
the spatial diagonal. By Proposition [HJ for such (z,z'), H n °(t, z, z') is phg conormal with 
infinite-order decay at \/i = 0, uniformly in the spatial variables; so V2(z)H no has the same 
properties. This proves the phg conormality part of the lemma, and the sup-norm statement 
is automatic since Qq is compact. 

It remins to consider Vi(z)H eZ (t, z, z')xi{z'). By the scaling property of the heat kernel, 

V x {z)W z {t,z,z')xi{J) = e- n V 1 {z)H z {t/e 2 ,z/e,z'/e)xM). 

We now use the coordinates (x,y,x',y r ). On eZ, V\{z) is a bounded linear combination of 
derivatives d r and d y , which are equivalent to rd r and d y since r G [9/8,5/4]. These lift to 
derivatives of the form xd x and d y on Z, so V\ lifts to a 6-differential operator on Z , which 
we call V\ as well. 

First we prove the phg conormality statement in the lemma. Fix an x and x' with x G 
(4/5,8/9) and x' G (16/17,2) and also fix y and y. Of course, e~ n is phg conormal on Q , 
so we only need to analyze 

V 1 H z (t/e 2 ,ex,y,ex',y'). 

We break this analysis into two cases: r = t/e 2 bounded above and r bounded below. 
First suppose that r < C . Then we are away from the face R of Q , so we use (a/t, e) as 
coordinates. We have: 

VxH z (t , ex, y, ex', y') = V 1 H z (r,ex,y,ex(x'/x),y / ). 

The spatial variables are fixed, so our variables are restricted to a plane in [0, C] T x Z x Z. 
Since x'/x is fixed and bounded away from 1, this plane is a p-submanifold of 

[Z 2 sc {x, y, x 1 , y') x [0, 1]^; { ^ = 0, (x, y) = (x', y')}]. 

The restriction H z {t, ex, y, ex' , y') is thus phg conormal on the plane; since x is bounded away 
from zero, (r, e) are good coordinates on the plane. We conclude that H z {t, ex, y, ex', y') is 
phg conormal in (r, e), and hence on Qq. Since V\ is a b-differential operator, all the same 
analysis holds for V\H Z . 

On the other hand, suppose that r > C^ 1 > 0. Then we use the coordinates (y/i,rj). 
Switching to the long-time heat kernel, we must consider 

V 1 F(ri,TiVix,y,7iVix> ',y'). (39) 

Our spatial variables are fixed, so we again restrict to a plane inside Z 2 SC ; this plane is a p- 
submanifold, and intersects bf, bf , and zf. Letting c = r/y/tx, we see that F(r], c, y, c(x'/x),y') 
is phg conormal on X 2 (i], c). Since we are concerned only with small t, c = r]\/tx < r//2, so 
F(r), c,y, c(x' /x),y') is phg conormal in the coordinates (i],c/r]) = (r],\/tx). However, since 
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x is fixed, F(r],r]\rtx,y,T]\/ r tx f ,y') is phg conormal in (r),\/i) and hence on Q . As before, 
V\ is a b-differential operator, so V\F is phg conormal on Qq as well. 

We must also check that the dependence on z and z' is smooth. However, this follows 
immediately from the fact that derivatives d r and d y lift to b-derivatives on Z, as we have 
already computed. 

Now we prove the supremum statement. Let N be any positive integer and consider the 
function 

rVTO) 8 M z (r, - e , ^) = r N e- n {W € {z')yV x H z {T, ex, y, ex', y'). (40) 

We will show that this function is globally bounded. By our support conditions, x G 
(4/5,8/9) and x' > 16/17. Therefore, the arguments of H z are bounded away from the 
'scattering' faces in both the short and long-time regimes, as well as rb and rbo for long 
time. Moreover, since t is bounded, r -1 / 2 > e > ex; so the arguments of H z are also 
bounded away from zf. Let x = ex' and e = ex. Then we have: 

r-^^TO)7rt $ V, x, y'). (41) 
We examine the behavior of (HIT) at each face of the blown-up (r, e, y, x, y') space from 
Figured First note that x is contained in (4/5,8/9) and hence is essentially constant. The 
function W e (z') lifts to be equal to e when x > 1 and 1/x' = x~ x {ejx) for all other x' (recall 
that x' > 16/17). It therefore has order 1 at the faces lb and lbo and zero at all other 
faces (it has order — 1 at rb and rb , but our support conditions avoid those faces). H z has 
infinite-order decay at bf, lb, and the t — face, and hence so does fl4"TT) . At bf , H z and 
hence V\H Z has order n, r~ N has order 2N ' , i~ 2N ~ n has order — 2N — n, and (W t (z')) a has 
order zero; the total order is thus zero. Finally, at rb , H z has order n — a, t~ n has order 
2N, i~ 2N ~ n has order -2N - N, and {W e {z')) a has order a, again giving a total order of 
zero. We conclude that (T4T]) is bounded at all boundary hypersurfaces which intersect the 
support of the arguments of H z , and hence that (14"T!) is globally bounded. Since N was 
arbitrary, this completes the proof. □ 

Now that we have analyzed E e , we form the Neumann series for K e and prove analogous 
short-time bounds for K e . Fix T > 0. 

Lemma 13. a) The expression E* k (t, z, z') is defined for each k, and the sum K e (t,z,z') = 
J2™ =1 (-l) k+1 E* k (t,z,z') converges. 

b ) For t <T and any MeN, there is a Cm independent of e such that 

sup \{W e {z')) a K e {t,z,z')\ < Ct M . 

z,z'GQ. e 

c) For each z, z' with r > 1/2, K e (t, z, z') is phg conormal on Q, with smooth dependence 
on z and z' in this region. Moreover, let the index set for E e (z,z',t) at e = be F = 
{{sj,Pj)}; then the index set for K e (z,z',t) at e = is F+, the set of all finite sums of 
elements of F. 

Note that F is bounded below by zero by Lemma [TTJ so F+ is an index set and is also 
bounded below by zero. 

Proof. \E e (z, z', t)\ <Cq uniformly for t in any bounded interval and z, z' G Q e ; therefore, all 
the integrals in the convolutions converge, and E* k (t, z, z') is defined. To see that the sum 
over k converges is slightly more involved. Since e is fixed for the moment, we suppress it 
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in the notation. All of the spatial integrands in the convolutions are zero when the variable 
of integration is not in the support of V\ or V2 (we call this the 'cutoff region'), so we can 
regard the integrals as being over a subset of Qq. Let V be the volume of Qq. Moreover, the 
weight function W e (z') is between 1/2 and 1 in the cutoff region. 

Now for any M and t in any bounded interval [0, T], Lemma ITU implies that there is a Cm 
so that \(W e (z')) a E(t, z, z')\ < Cm^ m ■ By integrating and using the weight function bound 
on the cutoff region, we get that 

\E* 2 (t,z,z')\ 



/ / E(t — s, z, z")E{s 1 z" , z) dz" ds 
Jo Jo.0 



< 



[ 2 a VC M (t - s) M C Q {W e {z')y a ds < C M t 
Jo 



,M (2 a VC t) k - 1 



We then use induction; we claim that \E* k (t, z, z')\ < C M (M\)t M ^M+k-iy (W £ (z'))- a for 
all k. Indeed, assuming the inductive hypothesis, we have: 



\E< k+1 \t, 



[ [ E(t- s, z, z")E* k (s, z", z) dz" ds 
Jo Jn 



< 
So: 



(2 a VC )C M (M\)s 







(M + k-l)\ 



(W e (z'))- a ds = C M (M\)t M{2aVCot)k 



\K(t,z,z')\<yC M (M\)t M{2aVCotri 



k=l 



(M + k-l)\ 



(W e (z'))- a < 



00 



/ CMt 
k=0 



(M + k)\ 

M ( 2 a VC t) k 
k\ 



(W e (z'))- a . 



(W e (z'))-° 



= C M t M e 2avCot (W e (z'))- a . 
This shows not only that the series converges, but also that for t < T, \K(t, z, z')\ < 
C M t M (W e (z'))- a for all M, where C M = C M e 2avc ° T ; since C M is independent of e, this 
completes the proof of parts a) and b). 

Now assume that z and z' both have r > 1/2. From Lemma ITTJ and the corollary following 
it, E e (t, z, z') is phg conormal on Q; we want to show that K £ (t, z, z') is also phg conormal 
on Q. Fix positive integers N and b. Since E e (t,z,z') is phg conormal, there exist I > 0, 
{si,Pi) G E for each % between 1 and I, and coefficients a{(t, z, z') for each i between 1 and / 

1 

E e (t, z, z') ~ <^ z ' z ') eSi ( io s e ) P! + R **& z ' z '^y ( 42 ) 



1=1 



Moreover, for each M G No, there exists C\,,m suc h that \di(t,z,z')\ < CMt M and Rn < 
Cb,Mt M e N , each together with up to b derivatives of the form (ed e ) or (td t ). 

When we take convolutions of E e with itself, we plug in ( FT2|) and work term-by-term. 
For each (r, q) G F + with 3?r < N, and each k, E* k (t, z, z') has a term of the form 
bk,r, q (t, z, z')e r (\oge) q ; bk, r ,q may be zero. We want to bound \bk,r t q(t, z, z')\. There are at 
most l k combinations of (si,Pi) with 1 < i < I which add up to (r, q), and these are the only 
combinations which contribute to bk^ r , q - For each combination, we can perform the same 
convolution analysis as before to bound the coefficient; we obtain that 

\b k , r , q (t,Z,z')\ < C b M M l k {V ^y - 
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The rest of the combinations yield terms of order e , as they either involve a remainder 
or their e exponents add up to more than N. We combine all of these terms into a new 
remainder term R^ t f.(t, z, z' , e). There are at most (I + l) k of these combinations. We again 
perform the convolution analysis to obtain that 

\R n>k {t, z, z>, e) | < C biM t M (l + l) ^^ ^" 1 ■ 

[k-l)\ 

Moreover, let V be any differential operator of order < b that is a combination of (ed € ) and 
(tdt) derivatives. We claim that 

\VR n . k (t,z,z',e)\ < C b>M t m (l + l) k k b i \^ b ' ot)k - . (43) 

(fc-l)! 

Indeed, each of the b derivatives in V may act on each of k different terms of the form 
cii(t, z, z')e Si (\oge) Pi or R N . If they act on R N , we still have the Cb t Mt M ^ N bound. On 
the other hand, if they act on aj(t, z, z')e Si (\oge) Pi , the ed e derivatives only improve the 
expansion, while the td t derivatives hit a,i(t,z,z'), and we can still apply the Cm^ M bound. 
This proves (l4"3l) . 

Finally, we sum from k = 1 to oo. By the Ratio Test, 



oo 



^ >r ,Q / , \bk,r,q(t } Z } Z ) | 
k=l 

converges and is 0(t M ) for any M. And for any differential operator V of order at most b 
as above, 



oc 



Rn = ^ ^ \R 



N,k\ 



k=l 



converges by the Ratio Test and is also 0(t e ). Therefore, we may write 

K e (t, z, z') = Yl &r,,e r (loge) 9 + R N . 

{{r,q)£T + ,$ir<N} 

Since b and N were arbitrary, K e (t, z, z') is therefore polyhomogeneous conormal on Q, with 
index set T+ at e = and infinite-order decay at \ft — 0. This completes the proof of 
Lemma [131 □ 



4.4. Construction of the heat kernel. We now analyze G e (t) * K e (t). If we can show 
that it is well-defined for each e, we will have shown that G e (t) — (G e (t) * K e (t)) is in fact 
the heat kernel H Q€ (t, z, z'). The next lemma proves what we need, and also gives a useful 
bound: 

Lemma 14. Fix T > 0. Let A t (t) be any family of operators on Q 6 whose Schwarz kernel 
A e (t,z,z') has z-support between r = 3/4 and r = 5/4, and suppose also that for each M 
there is a Cm such that for all e < 1/2, all z, and all t < T, 

\(W € (z')) a A £ (t,z,z')\<C M t M . (44) 

Then (G e * A e )(t, z, z') and (A e * G e )(t, z, z') are both well-defined for each e. Moreover, for 
each M , there is a Cm such that for all e < 1/2, all z and z' in Q e with z' between r = 3/4 
and r = 5/4, and all t <T, 

\(A e *G e )(t,z,z')\<C M t M . (45) 
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In particular, this lemma may be applied with A e = E e or A t = K e . 

Proof. First we consider A e * G t ; G e (t) has two terms, which we consider separately. The 
term coming from Q gives 

f I A(s, z, z")Uz")H n °(t - s, z", z') X2 {z') dz" ds. (46) 
Jo Jn 

Consider 

h 2 (s,z')= f Uz")H no (s,z",z') X2 (z')dz". (47) 
Jn 

By the definition of the heat equation, l%2(s,z) is precisely X2{z') times the solution of the 
heat equation on Q with initial data X2( z ")- That solution approaches X2{z') as s goes to 
zero and hence is uniformly bounded, for s < T, by a constant. So h 2 (s,z') is uniformly 
bounded for s < T. Combining this with (|44l) yields both the well-definedness and the bound 
for this piece (since z" is outside r = 1/2, W e {z") is between 1/2 and 3/2). 
On the other hand, the term coming from tZ is, substituting s = t — s, 

[ [ A{t- s, z, z")xi{z")H eZ {S, z", z')xi{z) dz" ds. (48) 
Jo Jn t 

We claim that 

f \(W e (z'')rU^W Z (^^,z')xi^)dz''\ (49) 

is bounded uniformly by some constant C for all z, all e < 1/2, and all s < T. Assuming 
the claim for the moment, the absolute value of (I48p is bounded, for any M, by 

C M t M [ C ds = CC M t, 
Jo 

which is more than enough. 

As for G e * A e , we only need well-definedness and not a uniform- in-e bound. So fix e and 
break G e into two pieces. The piece may be analyzed as above, switching s for t — s and 
X2 for X2- In the eZ piece, we get precisely 

f [ Xi(z)H £Z {s,z,z")x 1 (z")A(t- s,z",z') dz" ds. (50) 

Jo JeZ 

Notice that 

hi(s,e)= f Xi(z)H eZ (s,z,z")xi^") (51) 

JeZ 

is Xi{ z ) times the solution of the heat equation on the fixed manifold eZ with initial data 
Xi{z"). As before, the limit as s goes to zero of this solution is precisely Xi( z )i an d hence 
the solution is bounded uniformly for short time; we again use the known bound on A t to 
complete the proof. 

It remains only to prove the claim by analyzing (|49l) . We do this separately for z" inside 
and outside r = 1/2, which correspond to off-diagonal and near-diagonal regimes respectively. 
First suppose that z" is inside r = 1/2; remember that z' is between r = 3/4 and r = 5/4, 
so H eZ (s, z", z') is an off-diagonal heat kernel on eZ. Therefore, we may proceed precisely as 
in the proof of Lemma [Til with (z', z") instead of (z, z'), to conclude that for each M there 
exists Cm such that 

\H* z (s,z",z')\<C M t M (W e (z")) a 
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for all z" inside r = 1/2, all z' between r = 3/4 and 5/4, all e < 1/2, and all s < T. Plugging 
this bound into (1491) gives a bound for (|49j) of 

/ X i(^)(^")) 2a . (52) 

Since a < n/2, it is easy to compute directly that (15"2"j) is bounded, uniformly in e, which 
completes the proof of the claim for z" inside r = 1/2. 

On the other hand, suppose that z" is outside r = 1/2; in this region, W e (z") is between 
1/2 and 3/2 and hence we may ignore the factor of W e (z") a . We rescale the heat kernel and 
switch to an integral over Z; let %3 be equal to xi outside r = 1/2 and inside r = 1/2. 
Letting a = s/e 2 , z = z"/e, z = z'/e, x = and x = we get 

1/ / X3^/x)H z (a,x,y",x,y')xi(e/x) x^ 1 ' 1 dx dy"\. (53) 

J N Je/2 

The support conditions imply that x is between e/2 and 3e/2, while x is between 4e/5 and 
4e/3. So the ratio x/x is between 1/3 and 2. Let \a be a smooth cutoff, equal to 1 between 
1/3 and 2 and outside (1/4,3). The expression f l5"3"|) is therefore bounded, independent of 

e, by 

Xi(x/x) \H z (a, x, y", x, y') \x~ n ~ l dx dy". (54) 

Consider the integrand in ( 154"]) . It is phg conormal on the space in Figure [7]) , but it is 
also supported near the diagonal; that is, away from lb, rb, lbo, and rbo- It is simple to 
verify that integration in x is a pushforward onto the face rb crossed with N y . We apply the 
pushforward theorem of Melrose to conclude that (154")) is phg conormal on this space; the 
absolute value is irrelevant, even if the heat kernel is not assumed to be positive, because we 
are only interested in the leading orders. It remains to compute those orders. 

- The two faces of the space in Figure [7] which are mapped into a = are the a = 
diagonal and the a = non-diagonal. On the non-diagonal face, the integrand has infinite 
order decay. However, the integrand has order — n at the a = diagonal as a function. Since 
we are integrating with respect to dz = d{z — z") instead of d((z — z")/ ^fa), we get an order 
shift of n when applying the pushforward theorem. The result is that when applying the 
pushforward theorem we get a leading order of at the t = face of rb. (This order shift 
can be understood by considering the heat kernel on a compact manifold - the integral of 
the heat kernel in one of the spatial variables converges to 1 even though the heat kernel 
itself blows up to order n/2 at the t = diagonal). 

- As for the 0<cr<oo,x = face, bf and sc are mapped into it. The integrand has 
infinite-order decay at bf, and order with respect to scattering half-densities at sc, which 
are the appropriate half-densities to use there; the pushforward theorem gives a leading order 
ofO. 

-Atcr = oo,x = 0, the faces in the preimage are bfo and lbo- At bfo, the integrand has 
leading order n with respect to scattering half-densities, hence leading order with respect 
to the appropriate b-half-densities; since the integrand vanishes to infinite order at lbo, the 
pushforward theorem again gives a leading order of 0. 




28 



DAVID A. SHER 



- The faces zf and rbo are mapped into a = oo. The integrand is supported away from rbo 
and has leading order at worst at zf, so the pushforward leading order is 0. 

We conclude that (154"]) has non-negative leading orders at each boundary hypersurface of 
rbxiV^. It is therefore uniformly bounded, which proves the claim and completes the proof 
of the lemma. □ 

We can apply the lemma with A t = K e or A e = E e , so we conclude that all convolutions 
involving G and either K or E are well-defined. This argument completes the proof of 
Proposition [7J 

4.5. Proof of Theorem [31 Now that we have shown Proposition [7J we consider the trace 
of the heat kernel, which may be written as a sum of three terms as in (]3"Tj) . We now go 
term-by-term and prove that each is phg conormal on Q Q . In all cases, we restrict so that t 
is less than a fixed time T. First we show: 

Lemma 15. TrG t (t) is phg conormal on Q . Its leading orders on Q are at worst —n at 
L, —n at F, and at R. 

Proof. Tr G 6 (t) has two terms. The trace of the f2 term is just 

/ X 2(z)H n °(t,z,z) dz, (55) 

which is independent of e. Since the support of Xii?) ls aw ay from the conic tip, the integrand 
has a phg conormal expansion as t — > with smooth coefficients in z by Proposition [HJ 
Therefore, ( 15 5 p is phg conormal as t — > and independent of e, hence phg conormal on Q 
and therefore on Q . The leading orders of (155]) on Q are — n at y/t = and at e = 0, so 
the leading orders on Q are — n at L and F, and at R, precisely as needed. 
For the eZ term, we rescale, then change variables. Writing r = t/e 2 , we get: 

/ Xi{zW z (t, z, z) dz = e~ n [ xi(z)H z (t/e 2 , z/e, z/e) dz 

JtZ JeZ 

= [ X i(ze)H z (r,z,z) dz. (56) 
J z 

We observe that this integral is precisely ( 1201) . Recall from section 3 and [ST| that ( 1561) is 
phg in (r, e) for r bounded above and on X 2 (e, t~ 1 / 2 ) for r bounded below, which translates 
into polyhomogeneity in (r], y/t) as in section 3. Therefore (|56l) is phg conormal on Qq. It 
remains only to analyze the leading orders. This may be done by using the pushforward 
theorem and our knowledge of the leading orders of the heat kernel on Z to track the leading 
orders through the proof of the polyhomogeneity of (1561) (the proof of Theorem 12 in [Slj ) . 
For r less than any fixed time T, we do not need to use the pushforward theorem; instead, 
Theorem [TUl implies H z (t, z, z) < Cr~ n , and hence that 

J Xi(ze)H z (r, z, z) dz < CT~ n Vol{Z n {\z\ < 2/e}) = CY""^™. (57) 

Therefore the leading orders of (|56|) at L and F are each at worst —n. 

3 A key point in all these computations is that the integrand fails to decay to infinite order at only one of 
the boundary hypersurfaces in the preimage of each hypersurface of rbxiV^. The result is that there are no 
extended unions in the pushforward theorem and hence no logarithmic blow-up. 
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Now consider r greater than any fixed time T; we need to analyze the leading orders at R. 
The portion of the integral (|56|) with \z\ < 1 is a function only of r (as the cutoff function 
is identically 1 in this region for e < 1/2), and hence has leading order at R. We therefore 
restrict attention to the \z\ > 1 part of (I56p . which, as in the proof in |Slj . may be written: 



Here we must use the pushforward theorem. As in [Slj . we know that the integrand in 
f )58|) is phg conormal on X^(x,i],e). Let /ti be the leading order of F(i],x,y,x,y) at r\ = 
- that is, the leading order of the heat kernel on Z at zf - and let k 2 be the leading order 
of F(r],x,y,x,y) at rj — x — 0, which corresponds to bf . Then the leading orders of the 
integrand in ( JB"8l) . with respect to — , are oo at the x = face and at the x = r\ = face, 
k 2 — n at the x = rj = e = face, at the e = face, — n at the x = e = face, and Ki 
at the 7] = and 7] = e = faces. Applying the pushforward theorem shows that (1581) is 
phg conormal on X|(?7,e), with leading orders of — n at the e = face, k 2 — nOni at the 
e = rj = face, and n± at the rj = face; here the operation U is the extended union, which 
is discussed in [Mil IMe2llGr] . 

We may simply apply this analysis with the known values of k 2 = n and K\ — 0, but then 
we obtain a leading order of log e at e = rj = from the extended union; that face corresponds 
to R. To avoid this difficulty, we apply the remainder of the alternative hypothesis. Since 
the zero-order coefficient of F at zf decays to order strictly greater than n at bfo, we may 
write F — Fx + F 2 , where Fi decays to order at zf and n + 8± for some 5± > at bf , and 
F 2 decays to order 5 2 > at zf and n at bfo- Applying the above analysis to F± and F 2 
separately, the extended union problem vanishes, and we are left with a leading order of 
at e = T] = 0. As we claimed, this face corresponds to R; since e/r] = \/i is bounded, (|58|) is 
phg in (i], e/r]) = (r], with leading orders at r] = and — n at \ft = 0. Therefore the 
leading order of f )56|) at R is 0. This completes the proof of Lemma [15j □ 

We now consider (E e * G e )(t, z, z'); note that it is zero for z outside the cutoff region - that 
is, the union of the bands where V\ and V 2 are supported. We will prove the following: 

Lemma 16. [Technical Lemma] For z and z' in the support of V\ and V 2 , (E t * G e )(t, z, z') 
is phg conormal on Qo, with smooth dependence on z and z' . 

The proof of this lemma is quite long, and so it is deferred until the next section; we 
assume it for the moment. Once we have proven the lemma, though, it is easy to refine it 
further: 

Corollary 17. For z and z' in the support of V\ and V 2 , (E e * G t )(t, z, z') is actually phg 
conormal on Q, with infinite- order decay at t = and leading order at worst at e = 0. 

Proof. We obtain this result as an immediate consequence of Lemma fT6l and Lemma HU The 
bounds in Lemma O imply that (E t * G e )(t, z, z') has infinite-order decay at L and F and 
leading order at worst at R. In turn, since (E e * G e )(t, z, z') has infinite-order decay at the 
blown- up face in Qo, a standard argument implies that it is phg conormal on the blown-down 
space Q; see Proposition A. 3 in jS]. □ 

Now consider (K e * (E e *G e ))(t, z, z'); it is well-defined by associativity and Lemma [141 By 
Corollary [T7] and Lemma [TBI both K e and E e *G e are phg conormal on Q, with infinite-order 
decay at t = and leading order at e = 0, for z and z' in the cutoff region. We may 




(58) 
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therefore apply the same argument as in the proof of Lemma [13] to conclude that K 6 *E t * G e 
has the same properties. Finally, take the trace of both E e * G e and K e * (E e * G e ); restricting 
to the diagonal and integrating over the compact cutoff region, we have 

Corollary 18. Tr(E e *G e )(t) andTr(K e *(E e *G e ))(t) are phg conormal onQ with infinite- 
order decay at \/t = and leading order at e = 0. 

We have now shown that each of the three terms in (|37|) is phg conormal on Q . Therefore, 
their sum TrH Qc (t) is phg conormal on Q . Moreover, TrH ne (t) is TrG U( (t) plus a term 
which is phg conormal on Q with infinite-order decay at t = 0. This completes the proof of 
Theorem [3j even if we only assume the weaker 'alternate hypothesis.' 

Appendix A. Proof of the Technical Lemma 

In this section, we complete the proof of Theorem [3] by proving Lemma [T6l We decompose 
the convolution (E £ * G e )(t, z, z') into four terms: 

f [ V x {z)H' z (t - s, z, z") X i(z")H eZ (s, z", z')xi(z') dz" ds 

JO JeZ 

+ [ [ V 2 (z)H n °(t-s,z,z")x2(z")H n °(s,z",z')x2(z') dz" ds 
Jo Jn 

+ f I V^H^it-s^^^ixi&Xz^H^^z",^^ dz" ds 
Jo Jn 

+ f I V r 2 (z)^(t- S ,z,z")(X2Xi)(^)^ Z (s,^,^)Xi(^) dz" ds. 
Jo Jn 

The spatial integrals are all over Q e , but the first integrand is zero when \z"\ > 3/2 and hence 
the integral may be written as an integral over eZ; the remaining integrands are zero when 
\z"\ < 1/2 and hence their integrals may be written as integrals over Qo. We will examine 
each integral in turn and prove that it is phg conormal on Qo, uniformly for z and z 1 in the 
cutoff region. These proofs depend heavily on geometric microlocal analysis and Melrose's 
pushforward theorem. Some notes on the analysis to follow: 

- As before, we may treat V\{z) and V 2 {z) as cutoff functions, rather than differential op- 
erators; they lift to b-derivatives and thus have no effect on the polyhomogenous expansions. 
We suppress these cutoff functions in the integrals, as they do not interact with the variables 
of integration. 

- The support conditions imply that d^ f (z, z") > 1/16 whenever the integrand is nonzero, 
so z and z" are uniformly separated. On the other hand, z' may be close to either z or z" . 

- Since we are only interested in polyhomogeneity, we suppress all polynomial factors of 
the variables of integration; they may affect the leading orders, but not polyhomogeneity. 
All of the overall integrals, as well as all of the interior spatial integrals, are well-defined; as 
long as we avoid integrating the kernel on the diagonal in time down to t = 0, we will not 
have to worry about the integrability condition in the pushforward theorem. 

- Let p = s/e 2 , a = y/s/(rj\/i), and a' = (cr) _1 = (r]y/i)/y/s. It is often convenient 
to integrate in one of these variables instead of in s; the change of variables only gives 
polynomial factors, which we ignore. 
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- Many of the integrals will have to be broken into further pieces; among other decompo- 
sitions, we often have to distinguish between the situations when r is bounded above and 
when it is bounded below. 

A.l. The first integral. For the first integral, we rescale the heat kernels, then substitute 
z = z"/e so that we can integrate over Z instead of eZ. We get 

e -n /" f H 2 ^, z/e, z) X i(ez)H Z (^, 2, *7<0 d ~ z ds - ( 59 ) 
Jo Jz e e 

Using x and y instead gives, where x = ex": 

e-"£^F z (^,ex,y,£,y / Oxi(e/£)^ Z (^,£,^ea; / , 2 / / )^ 1 ^ dy" ds. (60) 

A. 1.1. Small-r case: First consider fl60l) where r < C for some C; in this case, since s < t, 
we use p = s/e 2 , and p < C also. Both heat kernels are thus in the short-time case. Rescaling 
the integral in s to an integral in p and suppressing polynomial factors, we obtain 

/ x({p<t}) H z (t - p,ex,y,x,y")xi(-)H z (p,x,y",ex',y') dx dy" dp. (61) 
Jo J z x 

The first heat kernel is at short time and away from the diagonal, so by Theorem [9J it is phg 
conormal on X%(ex, x) x [0, C]( r _ p ) x N y x N y »\ xi( e /%) is a l so P n g conormal on X 2 (ex, x), so 
its pullback is phg conormal on the same space. Since x and x' vary over subsets of [1/2, 2], 
the first heat kernel and the cutoff are also phg conormal on X 2 (ex', x) x [0, C\i T - p \ xN y x N y " 
with parametric dependence on x. The second heat kernel is, again by Theorem [9j phg 
conormal on [Z 2 c (x , y" , ex' , y') x [0,C] P ; {^fp — 0,x — ex',y" = y'}}. Therefore, the pullback 
theorem implies that the integrand in [61] is phg conormal on 

[Z 2 sc (x, y", ex', y') x [0, VC] ^ { y/p = 0, ~x = ex', y" = y'} x [0, C] (r _ p) , (62) 

with parametric dependence on (x, y). Integration in x, and y" is a 6-fibration from ( )62|) onto 
[0, l] ex ' x N y t x [0, l}^p x [0, 1] T __ P (see propositions A. 11 and A. 12 in jS], originally from |MS] 
and |Me2j ). Therefore, we may apply the pushforward theorem; the spatial integral in (!6T|) 
is phg conormal on [0, l] ex > x N y > x [0, C]^p X [0, C]( r _ p ) with (x, y) as parameters, and hence 
phg conormal on [0, 1] £ x [0,C] p x [0,C]( r _ p ) with (x,y,x',y') as parameters. 

All dependence on (x,y,x',y') is smooth, so we now suppress it. It remains to do the p 
integral. Notice that the map from [0, C] T x [0, 1] P / T to [0, C] p x [0, C] T _ P given by (x\, X2) — > 
(xiX2,x\(l — 22)) is a b-map. Therefore, the spatial integral above is phg conormal in 
[0, l] e x [0, C] T x [0, l] p / r . Change variables to an integral in pjr and integrate from to 1; 
this integration is a b-fibration onto [0, l] e x [0, C] T , so ( IBTl) is phg conormal in (e, r) and 
hence on Q . 

A. 1.2. Large-r case, introduction: Assume that r > C for some large C, and use (r),t) in 
place of (r, e). We break up the s integral into four pieces: I, from to e 2 = rft; II, from rft 
to t/2, IV, from t/2 to t — e 2 , and III, from t — e 2 to t. Each is a different regime. 

In regions I and II, the time argument of the first heat kernel is greater than t/2 > C/2, so 

we write it as F(^=, rjy/ix, y, x, y"), which equals F{V\I x- s /v ^V 7 ^, y, x, y"). We introduce 
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a placeholder variable a = r/y/t. The total integral in regions I and II is, again suppressing 
polynomials: 

--,ax,y,x,y")Xi(a/x)H z ( — ,x,y",ax',y') dx dy" ds. 
(s/t) rjH 

The support conditions guarantee that the arguments of F are away from sc. Therefore, 
F is phg conormal on X^(r/^j^,a,x) with smooth parametric dependence on (x,y,y"). 

However, ■sj it/ (t — s) is a smooth and bounded function of s/t, since s < t/2 in regions 
I and II. So: in regions I and II, F is actually phg conormal on X^(r],a, x), with smooth 
parametric dependence on (x,y,y") and also on s/t, which is between and 1/2. 

A. 1.3. Region I:. In region I, we have 



J J F (V]J 1 _ ^ t y ax,y,x,y'')xi(a/x)H z (-^-,x,y'',ax',y') dx dy" ds. (63) 

The temporal integral is from s = to s = rft, so the second term is a short-time heat 
kernel. We break the integral into two pieces, called I-a and I-b, with x < 2r]\/t and 
x > 2r]y/i respectively. 

Consider I-a; here x/(r)\rt) is between 1/2 and 2. The first heat kernel F is phg conormal on 
X^(rj, a, x), as is the cutoff xi {a fx). Since x/a is between 1/2 and 2, F is in fact phg conormal 
on X 2 (r], a) with parametric dependence on {x/iyjyi)) and s/t as well as on (x, y, y"). Using 
a = y/s/(r)y/t), the second term is phg conormal on the space 

[Z 2 c (x, y", V Vy, y') x [0, 1] CT ; {a = 0, x = ex 1 , y" = y'}]. 

Combining all the terms, the integrand is phg conormal on: 

Xl{a,rj) x ([Z 2 sc (£,y",vVtx',y') x [0, {a = 0, x = ex', y" = y'}}), (64) 

with parametric dependence on s/t = (rja) 2 as well as (x,y,y"). The spatial integral is in x 
and y", and the pushforward map is a b-fibration from fl64l) onto X 2 (a, r/) x [0, 1]-^./ x [0, l] a ; 
its fibers are transverse to the cutoff at x — 2r]\/t. By the pushforward theorem, the part 
of ( 1631) with x < 2i]\/t is phg conormal on X 2 (a,i]) x [0,1]^ x [0, 1] CT with parametric 
dependence on (x,y,x',y') and also on rj 2 a 2 . 

Now we consider the temporal integral; it goes from to 1 in a. We first need to restrict to 
a = r/y/t in the space X 2 (a, rj) x [0, x [0, 1] CT . X$(a, rj, rjy/i) x [0, l] a is a blow-up of this 
space, and restricting to a = r]\/i gives a function that is phg conormal on X 2 (i], r]\/i) x [0, l] ff 
with parametric dependence on (rja) 2 . However, this parametric dependence amounts to 
multiplying each coefficient by a smooth function of (rja) 2 ; we conclude that the spatial 
integral in (163]) is phg conormal on X 2 (r/, i]\/i) x [0, 1] CT . Integration in a is then a b-fibration 
onto X 2 (rj,r]yt). So the part of (1B"3"]) with x < 2r\yt is phg conormal on X 2 (r],r]\/t). 

So far we have treated as a formal variable, but now we remember that rjs/i/r) = \ft 
is bounded. Therefore we have an expansion in [r\, [r]\ft)/r\) = (rj, y/t), and the part of (163|) 
with x < 2r\\ft is phg conormal in [r\, \ft) and hence on Q . 

Now consider integral I-b. This is easier, as the second term is away from the diagonal. 
The first term is phg conormal on X^(r],a,x) with parametric dependence on (x,y,y") and 
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i/s/y/i = arj, and the second term is now phg conormal on X%(x,r)\rb) x [0, 1] CT with para- 
metric dependence on (x',y',y"). The cutoff functions are also phg conormal there. We 
immediately restrict to a = r/\/i and see that the integrand is phg conormal on 

Xjj( V ,riVi,x) x [0,1], 

with parametric dependence on (x, y, x', y', y") and rja. As before, the parametric depen- 
dence on (rja) does not affect the polyhomogeneity. Integration in x and y" is a b-fibration 
onto Xl(j],T)y/i) x [0, 1] CT with fibers transverse to the cutoff x = 2r)y/t, so we can use the 
pushforward theorem; the analysis then proceeds exactly as for integral I-a. Therefore ( 16"3l) 
is phg conormal on Qq. 

A. 1.4. Region II:. For region II, the third argument of the heat kernel is s/(r] 2 t), which is 
now larger than 1. So both heat kernels are long-time. We must analyze: 

/" <//2 f I i \ft 

I / F (W 7 > rjVix, y, x, y")xi{r]\/t/x)F{r]—, x, y", rjVtx', y') dx dy" ds. (65) 

Jrft Jz V t-s Vs 

As before, break the x integral into two pieces: Il-a, where x < 2r)y/i, and Il-b, where 
x > 2i]\/i. In each pieces, switch to an integral in o' = (r]\/i/^/~s). In both pieces, the first 
kernel is phg conormal on X^(r], r)\rt, x), but with parametric dependence onrj/a'] hence the 
first kernel is phg conormal on X£(r], r]\/i, x, a'). It has parametric dependence on (x, y, y"). 

Consider region Il-a. As in the analysis of I-a, the restriction to x < 2rj\/i means that 
x/(r)y/t) is between 1/2 and 2. We again introduce the placeholder a, but this time we also 
introduce a placeholder b for a'. Since x/(r)y/t) is bounded, the first term is phg conormal 
on Xjf(a, T], b) with parametric dependence on x/ (rj\/i) as well as (x, y, y"). The cutoff is phg 
conormal on X%(x,rj\/t). The second kernel is phg conormal on Z 2 , sc (er', x, y", T]y/tx\ y'). 
Therefore, the whole integrand is phg conormal on 

X*(a, V ,b) x Z 2 a , sc (a',x,y",r ] Vix',y'), (66) 

with parametric dependence on (x,y). Integration in x and y" is a pushforward by a b- 
fibration onto Xj*(a, 77, b) x X 2 (a', r]yt). So the spatial integral in ( 1631) over x < 2r)y/t is phg 
conormal on Xj}(a,r},b) x X^io' ,rj\/i), with parametric dependence on (x,y,x',y'). By the 
pullback theorem and the theory of b-stretched products in |MSj (also see the remarks after 
proposition A. 10 in [S]), it is phg conormal on X^(a,r],b,a',r]\/t). 

We need to resolve the placeholders; restricting to a = r/\/i and then b = a', the result 
is phg conormal on X^(r),a',r)y/i). We now integrate in a' from rjy/2 to 1; this integral is 
a pushforward by a b-fibration onto X%(r),r) So the part of ( 1631) with x < 2r]\/t is phg 
conormal on X 2 (r/,r]\/t), and thus phg conormal on Q as in the analysis of integral I-a. 

As for Il-b, the first kernel is phg conormal on X*(x, T], rjyt, a 1 ). The second kernel is now 
supported away from the scattering diagonal, and hence is phg conormal on X^(a', x, rjyt); 
and the cutoff is phg conormal on X%(x, r/y/t). By the pullback theorem, the integrand is phg 
conormal on X£(x, i], i]\/t, a') with parametric dependence in (x, y, x', y', y"). Integration in x 
and a', as before, is pushforward by a b-fibration onto X 2 (r], rjyt), and the analysis proceeds 
as in Il-a. We conclude that (1631) is phg conormal on Qq with smooth dependence on z and 



z'. 
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A. 1.5. Region III:. In both region III and region IV, the second term is a long-time heat 
kernel. Make the substitution s = t — s; our integral becomes (modulo polynomial factors) 



t/2 



HZ (^,vVtx,y,x,y")xi(vVi/x)F(r}\ -^-,x,y",r]Vix',y') dx dy" ds. (67) 
rjH V t — s 



First consider region III, which is the integral (|67|) from s = to s = rft. Change variables 
to a — \Z& j (r]\rt); it is integrated from to 1. H z is a short-time heat kernel away from 
the diagonal and hence is phg conormal on X^(rj\/i, x) x [0, 1] CT , with parametric dependence 
on (x,y,y"); Xi(.V\/t/%) is phg conormal on the same space. Therefore, the first heat kernel 
and cutoff are also phg conormal on X^(i]\/tx',x) x [0, 1] CT , with parametric dependence on 
(x,x',y,y"). 

As with the integrals in region I, F is phg conormal on Z^ sc (r], x,y" ,r]\/tx' ,y'), with 
parametric dependence on a/I/a/t^ = r\a between and 1/2 (the parameter isn't involved in 
the spatial coordinates on Z^ sc ). Let c be a placeholder for rja. By the pullback theorem, 
the integrand of ( 1671) in region III is phg conormal on 

Zl sc ( V ,x,y", V Vix',y') x [0,1]* (68) 

with parametric dependence on (x,y,x',c). The spatial integral is a pushforward by a b- 
fibration onto X^(r], r]\/tx') x [0, 1]*, which is the same as X£(r), ri\fi) x [0, 1]*, with parametric 
dependence on (x, y, x,' ,y', c). So in region III, the spatial integral in ( 1ST1) is phg conormal 
on X^(r],r]\/i) x [0, 1} 5 with parametric dependence on c = (na) as well as (x,y,x',y'). We 
now follow the analysis of integral I-a, with o in place of cr, to conclude that the part of flBTl) 
corresponding to region III is phg conormal on Q . 

A. 1.6. Region IV:. Finally, we have region IV, which corresponds to s — rft to s = t/2. 
Both heat kernels are long-time, and we need to analyze: 

/ F (V^, rjVtx, y, x, y")xi(j)\ft/x)F{r)J - — -, x, y" , rjVtx', y') dx dy" ds. (69) 
V H Jz Vs V t-s 



As before, change variables to a' = r\\ftl\f%\ then a' goes from rj/y/2 to 1. The first term 
is an off-diagonal long-time heat kernel, and is thus phg conormal on X^(a' ,r]\/t, x) with 
parametric dependence on (x, y, y"); so is the cutoff function Xi( ? 7V / ^/5)- The second term is 
again phg conormal on Z^ sc (r], x, y", r]\/ix', y') with parametric dependence on a/I/ \ft = r]/a' 

smooth down to Vs/Vt = 0. Introduce three placeholders a = r]\/tx', b = x, d = r\. Then 
the integrand of flBT?]) is phg conormal on 

X*(a', a, b, d) x Zl sc ( V , x, y", V V~tx' , y') (70) 

with parametric depndence on (x, y, x', y"). 

All the integrals are absolutely convergent since a' does not go down to zero. So we 
can switch the integrals and integrate in o' first. The a' integral is phg conormal on 
X^(a,b,d) x Z^ sc (r],x,y",r]\/tx',y'). Since a = r]y/tx', b = x, and d = 77, we would like 
to say that the (a,b,d) part is phg conormal on X^(r]\/ix', x, 77) and hence on its blowup 
Z^ sc {ri ) x 1 y" ) rj\/ix\y'). We could then conclude that the whole a' integral is phg conormal 
011 z lsc(v,S:,y",vVix',y'). 
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To justify this assertion, note that the a' integral has expansions at all boundary faces 
and corners of Z^ sc (r],x,y",r]y/tx' ,y'), with coefficients in Ap hg (X^(a, b, d)) for some fixed 
index family S on Xf. It may therefore be written as a sum of terms of the form Wjfj, where 
Ui is phg conormal on Z^ sc (i], x, y", r/y/ix', y') with index family Ti approaching infinity, and 
Vi is phg conormal on X^(a,b,d) with fixed index family S. Replacing the placeholders, 
Vi is phg conormal on X^(rj,x,r]y/tx') with fixed index family, and hence is phg conormal 
on the pullback Z 2 c (r),x,y'',r)\/tx',y') with fixed index family £K Therefore u{Vi is phg 
conormal on Z 2 sc (r),x,y'',r)y/tx',y') with index family £" + J P i, which approaches infinity 
as i increases. The sum over i is therefore phg conormal on Z 2 c (rj,x,y",r)y/tx',y'), with 
parametric dependence on (x,y,x r ). 

In conclusion, ( l69p is the integral over x and y" of a function which is phg conormal on 
Z^ sc (r], x, y", rjy/tx', y'). By the pushforward theorem, (1551) is phg conormal on X 2 (rjy/tx', rf) 
and hence on X 2 (r]y/t,rj) with smooth dependence on (x,y,x',y'). By the same analysis as 
in region I, (1691) is thus phg conormal on Q . This completes the analysis of region IV, and 
hence of the first integral. 

A. 2. The second integral: The second integral is 

f f V 2 (z)H n °(t-s,z J z")x2(z")H n °(s,z",z')x2(z') dz" ds. 
Jo Jn a 

This integral is independent of e. Moreover, by the same arguments as in the proof of Lemma 
[HI it decays to infinite order in t at t — 0, uniformly in the spatial variables as they range 
over the cutoff region. Therefore, this integral is phg conormal on Q and hence on Qq, 
uniformly in z and z'. 

A. 3. The third integral: The third integral, up to polynomials in e and suppressing V\(z) 
and X2(z'), is: 

f I HZ^Ai-Kxim^H^z",*') dz" ds. (71) 
Jo Jn a e e e 

Note that now z" is supported between r = 1/2 and r = 2. We again break it into the 
small-r and large-r cases. 

A. 3.1. Small-r case: For small r, we let p = s/e 2 and then let /x = 1 — p/r. Suppressing 
polynomials, we get: 

f H z (fi,-,-)(xiX2)(z")H Q °(^T(l - n),z",z') dz" dp. (72) 
Jn e e 

The first heat kernel is off-diagonal, and hence is phg conormal on [0, 1] M x [0, l/2) e with 
parametric dependence on z and z". Let a be a placeholder for e 2 r(l — p); the second term 
is on-diagonal but away from the conic tip, and hence is phg conormal in (z", z', y/a) with a 
blowup at {y/a = 0, z" = z'}. The whole integrand is phg conormal on 

[[0, 1], x [0, l/2) e x [0, 1]^ x (z\ z")- {v^ = 0,z' = z'% (73) 

with parametric dependence on z. Integration in z" is a b-fibration onto [0, 1] M x [0, l/2) e x 
[0,1] ^5, so the spatial integral in ( 172]) is phg conormal on [0, 1] M x [0, l/2) e x [0,1]^ and 
hence on [0, 1] M x [0, l/2) 6 x [0, l] a , with parametric dependence on (z, z'). 
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The map from [0, 1]„ x [0, l/2) e x [0, VC]^ to [0, 1] M x [0, l/2) e x [0, l] a given by (//, e, r) -> 
(//, e, ey / r(l — //)) is a b-map, so by the pullback theorem, the spatial integral in (172|) is 
phg conormal on [0, 1]^ x [0, l/2) e x [0, \fC\^f. Integration in /i is then a b-fibration onto 
[0, l/2) e x [0, v^C*]^, and hence f!72|) is phg conormal in (e, -y/r) and therefore on Qo- 

A. 3. 2. Large-r case: Here we use (77, Vt), and switch s = t — s. Break the integral into two 
integrals, one from s = to s = rft and the other from rft to t. 

For the first integral, let a = \f&/(r]\/t), and the integral becomes (up to polynomials): 

f [ H z {a\ 4=, ^)( Xl x 2 )(z")H Q o(t(l - a 2 r, 2 ), *) dz" da. (74) 

The first term is off-diagonal and hence is phg conormal in (a,r]\/t), and therefore phg 
conormal in (a, r), with parametric dependence on z' and z" . As in the small-r case, let 
a be a placeholder for t(l — a 2 r] 2 ), and then the integrand is phg conormal on 

[[0,l},x[0,C} v x[0,VT} v - t x[0,l] VE x(z',z");{^ = 0,z' = z"}} (75) 

with parametric dependence on z. As before, the spatial integral in (IT4"|) is phg conormal on 
[0, 1]^ x [0, C] v x [0,T] t x [0, l] a , with parametric dependence on (z,z'). Therefore it is also 
phg conormal on [0, l] s X [0, 1/2],, x X 2 (t, a). 

Since 1 — o 2 rf is bounded away from zero for small 77, the set {a = t(l — a 2 i] 2 )} is a 
p-submanifold of [0,1]^ x [0,1/2],, x X 2 (t,a). By the usual restriction theorems for phg 
conormal functions (see |Me2] or proposition A. 8 in [S]), the spatial integral in (1741 is phg 
conormal on [0, \] 5 x [0, 1/2],, x [0, Integration in a is a b-fibration onto (rj, space, 
and therefore ( 1741) is phg conormal in (77, \fi) and hence on Qq. 

In the second integral, the first heat kernel is long-time, so we switch to F. Switch to 
o' = rjy/i/yfi, and the integral becomes, modulo polynomials: 

f X({V < a'}) [ F(a', V Vtx,y, V Vix",y")H^(t(l - ^) 2 )), z ", z >) dz" da'. (76) 
Jo Jn ° 

The first heat kernel is off-diagonal, hence phg conormal on X 2 (a', rjyt), with (z, z") param- 
eters. Let a be a placeholder for t(l — (r//cr') 2 ). Then exactly as in the analysis of (I74|) . the 
spatial integral in ( !76|l is phg conormal on X 2 (a' ,r]\/t) x [0, 1] , with (z,z r ) parameters. 

Let X be the portion of X%(a', rj, r]\/i) where a' > r\ and [r]\/i)/r] < \/T, with a boundary 
face at a' = r\. Then X is precisely [0, VT]^ x [0, l] v/9 , x [0, 1]^. The map from X to 
X 2 (a',r]^t) x [0,l] a given by (y/i,rj/a',a') -> (a, (y/i(ri/a')a'), y/i(l - {ri/a') 2 )) is a b-map. 
Moreover, the cutoff x({&' > v}) is a l so P n S conormal on X. Therefore the whole a' integrand 
is phg conormal on X and hence on Xjj(a', rj, r)\/t). We now follow the analysis of region 
Il-a of the first integral to conclude that ( ITS"]) is phg conormal on Qq. This completes the 
analysis of the third integral. 

A. 4. The fourth integral: The fourth integral, suppressing Vi{z) and Xi{ z ') and again 
discarding polynomial factors, is: 

f I H^(t-s,z,z")(X2Xi)(z")H z (^,-X) dz" ds. 
Jo J no e 2 e e 
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The H n ° term is away from the diagonal; it simply decays to infinite order in (t—s), uniformly 
in z and z" . The calculations are therefore a much-simplified version of those in the analysis 
of the first integral. 

A. 4.1. Small-r case: First assume that r is bounded; let p = s/e 2 , and we get, suppressing 
all cutoffs and polynomials: 

f f H n °(e 2 (r- p),z,z")H z (p,x,y'\ex',y') di dy" dp. (77) 
Jo Jz 

. The first heat kernel is phg conormal on [0, l] e x [0, l] T _p, with parametric dependence on x, 
y, x, and y" . Therefore, the whole spatial integrand is phg conormal on [Z 2 c (x , y" , ex' , y') x 
[0, l] p ; {y/p = 0, z" = z'}} x [0, l] T -p with parametric dependence on z. By the pushforward 
theorem, the spatial integral in (1771) is phg conormal on [0, l] p x [0, l] ex i x [0, l] T - p . As with the 
short-time piece of the first integral, it is actually phg conormal on [0, 1] T x [0, l] p / T x [0, l] e ; 
by the same analysis, (177)) is phg conormal on Q . 

A. 4. 2. Large-r case: Now assume that r > 2. Break the s-integral into two pieces: s < rft 
and s > rft. For the s < rft portion, we have: 

I [ H n °{t{l -a 2 r ] 2 )),z,z")H z {a 2 ,x,y",r ] Vix',y') dx dy" da. (78) 
Jo Jz 

Since rj < 1/2 and a < 1, the first kernel is a smooth function of (t, a, rj), and hence smooth 
on X 2 (i]\/i,r]) x [0, l] ff (remember that \ft = (rjy/i/rf) is bounded). As in part I-a, introduce 
the placeholder a for rjy/t. The second kernel is phg conormal on the same space as in the 
small-r case above, with r]\/i replacing e. Therefore, the whole spatial integrand is phg 
conormal on 

[Z 2 C (5;, y", V Vix', y) x [0, l] ff ; {a = 0, *" = z'}} x X 2 (a, V ) (79) 

with parametric dependence on z. Integration in (x, y") is pushforward by a b-fibration onto 
X 2 (a, rj) x [0, 1 ]—/£,./ x [0, l) a . From here, we follow precisely the analysis in I-a to show that 
(178)) is phg conormal on Q . 

Finally, consider the integral from rft to t. Using a' instead, we have 

/ / H Qa {t{l - {v/cr') 2 )), z, z")F(a', r/Vix" , y", rjVix', y') dx" dy" da'. (80) 

Jrj Jz 

We proceed analogously to the calculation in region Il-a; the first term is a smooth function 
of \/i and r//a', so by the pullback theorem, it is phg conormal on the region of X^rjy/i, r], a') 
where (r]\/i/r] = y/t) is bounded and rj/a' is bounded, and hence on Xjj(r)\/i, 77, a'). As in 
Il-a, introduce placeholders a for r/y/t and b for a'; the analysis is then identical to Il-a, and 
we conclude that (180 j) is phg conormal on Q . This completes the analysis of the fourth 
integral, and with it, the proof of the technical lemma. 
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